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Abstract 

We introduce and study the properties of the "color-straight" four-quark oper- 
ators containing heavy and hght quark fields. They are of the form {hTi,h){qTqq) 
where both brackets are color singlets. Their expectation values include the bulk of 
the nonfactorizable contributions to the nonleptonic decay widths of heavy hadrons. 
The expectation values of the color-straight operators in the heavy hadrons are 
related to the momentum integrals of the elastic light-quark formfactors of the re- 
spective heavy hadron. We calculate the asymptotic behavior of the light-current 
formfactors of heavy hadrons and show that the actual decrease is l/(g^)2 rather 
than The two-loop hybrid anomalous dimensions of the four-quark operators 

and their mixing (absent in the first loop) are obtained. Using plausible models 
for the elastic formfactors, we estimate the expectation values of the color-straight 
operators in the heavy mesons and baryons. Improved estimates will be possible in 
the future with new data on the radiative decays of heavy hadrons. We give the 
Wilson coefficients of the four-fermion operators in the l/m^ expansion of the inclu- 
sive widths and discuss the numerical predictions. Estimates of the nonfactorizable 
expectation values are given. 
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1 Introduction 

The heavy quark expansion proved to be useful in describing decay properties of 
beauty hadrons. At the level of nonperturbative effects a number of local heavy 
quark operators of increasing dimension appears whose expectation values in the 
heavy flavor hadrons determine the importance of preasymptotic effects. The first 
nontrivial operators, chromomagnetic Oq = b^cr^j_yG^uh and the kinetic operator 
0,r = h{iD)% have D = 5. The expectation value of Oq is known directly from 
the masses of beauty hadrons. The expectation value /i^ of Ot, is not yet known 
definitely, although a certain progress has been achieved over the last few years in 
evaluating it for B mesons. 

More operators appear at D = 6, in particular, four-fermion operators hVhqV'q 
where q are light quarks and F, V denote various Lorentz and color structures. In the 
inclusive widths of heavy hadrons such expectation values govern 1 / m\ corrections. 
Their effect is still significant, especially due to specific accidental suppression of the 
impact of the leading D = 5 operators. 



1 



Unfortunately, the expectation values of the four-fermion operators up to now 
remain rather uncertain. Since the mid-80s |]T|, the vacuum factorization approxi- 
mation has been used to estimate the mesonic matrix elements which then appear 
proportional to Such factorizable terms are absent in baryons, and a number of 
simple constituent quark model estimates have been employed. The validity of the 
assumptions implemented in such analyses is not clear, however. As a result, the 
expectation values where the factorizable contributions are absent or suppressed, 
remain uncertain. 

On the other hand, the problem of a more reliable evaluation of the relevant 
four-fermion expectation values recently attracted a renewed attention since the 
lifetime ratios of the different beauty hadrons have been accurately measured. While 
data and predictions of the meson lifetimes are non-trivially consistent, the small 
experimental ratio taj^/t^^ = 0.78 ± 0.07 [3, if taken literally, seems to be in a 
conflict with the expectations based on the l/mt expansion. 

In the framework of nonrelativistic quark description the four-fermion expec- 
tation values are all expressed via the wavefunction density at origin |\E'(0)p (for 
mesons) or the diquark density / d^y\'^{0,y)\'^ (for baryons). All expectation values 
differ then by only simple color and spin factors For example, in B mesons one 
has 

^ ■{B-m{uu)\B-) = |M/(0)p, 



2Mb 

1 



{B-\{M-f,u){m-f,b)\B-) = N,\^{0)\' (1) 



2Mb 

(color indices are contracted inside each bracket), and for baryons 

1 



2Ma 



■{A,\ibb){uu)\A,) = d'y\'^iO,y)\ 



■{M\{bu)iub)\A,) = \ Jd'y\^{0,y)\\ (2) 



2Mb 
etc. 

In actual QCD this simple picture does not hold, and the naive Quantum Me- 
chanics (QM) relations between different expectation values are generally violated. 
Moreover, the notion itself of the nonrelativistic wavefunction used in the potential 
description, becomes ambiguous. Even in the perturbative domain the expectation 
values become scale-dependent, and the renormalization is in general different for 
different operators. This manifestly goes beyond the potential description, even 
extended for the price of introducing various light-quark spin wavefunctions in an 
attempt to account for the relativistic bispinor nature of the light quark fields. 

In this paper we note that there exists nevertheless a natural generalization of 
the notion of the wavefunction density, in particular at origin (the origin is defined 
as the position of the heavy quark). It is associated with the expectation values 
of those four-fermion operators for which the bTb bracket is a color singlet. The 
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color flow for such operators is not disturbed, and we call them "color-straight" 
operators. Their expectation values in the heavy quark limit rrih oo are related 
to the observable transition amplitudes. This fact suggests that they are better 
candidates for the operator basis used to parametrize hadronic expectation values 
in various applications. Moreover, they are more suitable also for applying general 
bounds of the type discussed in 0]. Such QM-type inequalities can be formulated 
more rigorously for these operators in full QCD. 

Knowledge of the light-quark current elastic formfactors of heavy hadrons would 
allow one to determine the color-straight expectation values. Unfortunately, they 
are practically unknown yet. Nevertheless, employing reasonable assumptions about 
their dependence allows more definite estimates of the expectation values. As the 
most conservative attitude, they can be viewed as educated dimensional analysis, 
with the added bonus of being free of ambiguities related to ad hoc powers of 27r 
inherent in various naive dimensional estimates. Such numerically significant un- 
certainties often cause controversy in the resulting expectations leading sometimes 
to rather surprising phenomenological conclusions. We also think that the derived 
relations can be used for an alternative, simple evaluation of the color-straight ex- 
pectation values in the lattice heavy quark simulations. 

2 Color-straight operators and light current am- 
plitudes 

Our main object of interest is the expectation values of the color-straight operators 
of the generic type 

hThb'q,Tgq^ (3) 

where r^, Tg are arbitrary matrices contracting Lorentz indices (Fg can be also a 
matrix in the light flavor space), and i, j are color indices. We will consider the heavy 
quark limit mf, oo assuming that the normalization point /i of the operators or 
currents is set much smaller than m^. In this case there are two nonvanishing types 
of operators transforming under rotations of the heavy quark spin as spin-singlet 
and spin-triplet, respectively. The corresponding Dirac structure on the heavy side 
is Fb = 1 and F;, = 775 = a : 

0._. = [bb) (gF,g) , O..,. = (gF,g) (4) 

All possible F{,-structures are reduced to these operators. In our discussion we always 
assume that the heavy quark is at rest, = (1,0), and denotes the velocity of 
the 6-hadron Hi,. Since in the heavy quark limit the b quark spin decouples, we start 
for simplicity from considering the spin-singlet operators O^.^. The straightforward 
generalization for 0^_^, will be formulated later. 

If a heavy meson were a two-body QM system where, additionally, the light 
quark is nonrelativistic as well, the expectation values of Oj measure the meson 
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wavefunction at origin, see the first of Eqs. (|l|), and likewise for otfier matrices Tq for 
which different spin wavefunctions ^(x) can enter. In the momentum representation 

(we use the normalization where / d^p/(27r)'^ |\E'(p)p = 1). 

On the other hand, in such a nonrelativistic system the Fourier transform of the 
light quark density distribution measures the elastic transition amplitude (formfac- 
tor) of the meson associated with the scattering on the light quark: 

^Q) = ^{B{q}\qq{0)\Bm = J d^f ^(f)^*(x)e "^^"^ . (6) 

The following relation then obviously holds: 

/ (0-^^^") = = ^{B\m{mm\B) . (7) 

Integrating the transition amplitude over all q yields the local four-fermion expecta- 
tion value we are interested in. Since we study a transition induced by scattering on 
the light quark, the scale of the transferred momentum is the typical bound-state 
momentum and is much smaller than m?,. 

In actual QCD the simple nonrelativistic picture does not apply. The light quark 
is certainly relativistic. Additionally, a two-body potential description (generally, 
any fixed-parton wavefunction) can only be approximately correct, with a priori 
unknown accuracy. 

It appears, however, that in spite of the fact that neither Eqs. (|]) nor (^ can 
be rigorously written in QCD, the final relation between the momentum integral of 
the (elastic) transition amplitudes and the color-straight expectation values holds 
exactly, up to corrections vanishing when mi, oo. It is not difficult to see, for 
example, that proceeding from a two-body nonrelativistic meson to a three-body 
nonrelativistic baryon does not modify the relation. We do not illustrate it here, 
and instead give a general field-theoretic proof. 

Let us start with the operator Og-s = hhqVq and consider the corresponding 
light quark current and its transition amplitude: 

Jr(x) = qTq{x)- -l-{H,(q)\J^(0)\H,iO)) = Ariq) • (8) 

The current does not need to be scalar; any particular component can even be 
considered separately. Likewise, the transition amplitude may not be a true scalar. 
The initial and final states may differ. The following relation holds: 



^j^-{HmbbqTqmHm = J Mo) ■ (9) 
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To prove this relation, we write explicitly the state Hb{q) with non-zero momen- 
tum as a result of the Lorentz boost from rest to the velocity v = q/Mjj : 



U 



L 



(10) 



where U 



is the corresponding Lorentz boost unitary operator. This op- 



erator is given by |^ 

U[L{v)] 



-in-Kd 



sinh^ 



n 



V 

\v\ 



the boost generators K can be expressed in terms of the symmetric energy-momentum 
tensor T^^: 

K' = j d^f (x'r°° - x°T°') (12) 



[x^ is fixed in Eq. ( [TI1| ) and can be put to zero). Since q does not scale with m^, 
we actually need to retain only the linear in v terms, which leads to simplifications. 
For example, the polarization degrees of freedom of (if any) do not change at the 
boost. 

The whole energy-momentum tensor consists of two parts: 

T,u = Tlf' + Tl^y = TX' + \b[l,m. + lum,- (zB).7.-(%7.]fc, 

(13) 

where Tj^^^^ is the usual QCD energy-momentum tensor including only light fields; 
it is free of the large parameter rub- In the heavy quark limit we need to retain only 
the part of T^i, which is proportional to m^: 



hh{x) + O ( 



(14) 



Here we have used the equations of motion for the h field. The anomalous terms are 
included in the last term (see, e.g., Sect. II). Therefore, we arrive at 



[HbmjvmHbm = {Ht 



The heavy quark limit leads to further simplifications: the number of heavy 
quarks becomes fixed and b itself becomes static. Then in the single-6 sector the 
following identity holds: 



,i / d''zf{z)bb(z) 



bb{z) . 



(16) 



Indeed, in the single-6 sector any product of the static b quark bilinears is very 
simple: 



{baib){zi)...{banb){zn) = 6^{zi - Zn)...S^{zn-i - Zn) bai...anb{zn) 



single b 



(17) 
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(cT/fc are arbitrary spin matrices). Using this, we obtain 



n=0 

oo -n 



J2 - r{z)bb{z) = / dh e'f^'Hb{z) . 

n=0 •' 

Taking f{z) = qzwe rewrite Eq. (|T5|) in the desired form: 

(^,(g)|Jr(0)|i/,(0)) = {mm J d'ze''^'bb{z)MO)\H,{0)). (19) 

Eq. (^) provides the discussed quantum field-theory generahzation of the notion 
of the hght-quark density qTq at arbitrary separation; one can define, for example, 

'"^"^"^^'^^ = / (0'^'''2ii^^^''^^''^''^^^°^'^^^ 

In what follows we are interested in the local heavy quark operators, that is 
when the light field operators enter at the same point as the b quark field. It is these 
operators that appear in the heavy quark expansion. Integrating Eq. (plOf ) over q we 
get 

{Momo)JrmHm = J mmrrnHm ■ (21) 

This is our master equation. We see that, in principle, it is even more general 
than was stated earlier: Jr(0) can be arbitrary gauge-invariant operator composed 
of the light fields, and not necessarily a light-quark bilinear. Besides, this relation 
holds not only for the truly forward transition matrix elements. The initial and final 
state hadrons can be different. Generally, they can even have different momenta; 
however, it must be assumed that these momenta are small compared to mi, - say, 
of the typical light hadron mass scale. Since this equation involves the integration 
over all transferred momenta, varying the relative momentum of the final and initial 
hadrons have no effect whatsoever, as it should be. 

Informative relations emerge, on the other hand, if we vary the heavy flavor state 
\Hf,) (or \Hf))) within the corresponding heavy-spin multiplet. Since the 6-quark spin 
decouples, this yields similar relations for the color-straight spin-triplet operators 
containing bab, that is, with the axial- vector 6-quark current. In particular, 

{Ht,mb^kbiO)JrmHbm = J -0^{SkHtmJrmHtm , (22) 

where 5/2 is the 6-quark spin operator. Formally one obtains this using, for example, 
the representation 

\SkHb) = J d'xbakbix)\Hb) 
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and applying relations ([T6|), (|17D generalized to include the 6-quark spin matrices. 
Alternatively, it follows merely from the heavy-spin symmetry relation between ma- 
trix elements of the operators 66Jr(x) and bcrj^bJrix) . 

It is worth to give a less rigorous but a transparent QM derivation of the master 
equation Eq. (^). Let us represent the expectation value of the color-straight 
operator bsb Jy{0) {s is either the unit or a spin matrix) by the sum over possible 
intermediate states: 

(^,M(0)Jr(0)|i^,) =E / {H,\bsbmn{q)){n{q)\.MO)\H,) . (23) 

The states \n{q)) are hadrons with a single b quark. In the effective theory the 
integral over momenta must converge at a hadronic scale which is much smaller 
than mfc. Then only the elastic transition (i.e., where and \n) belong to the same 
hyperfine multiplet differing, at most, by the heavy quark spin alignment if s is not 
a unit matrix) survive in the sum: all excited transition amplitudes generated by the 
heavy quark current 6s6(0) are either proportional to 1/m, or to velocity v ~ q/rrib 
of the heavy hadron state \n{q)) Moreover, since — ^ the elastic amplitude is 
unity up to corrections ~ q'^/ml we neglect. Thus, Eq. ( plf ) is reproduced. 





Figure 1: Diagrams for the renormalization of the four-fermion operators. 

Let us illustrate the validity of relation ( ^I]) diagrammatically, in respect to the 
perturbative corrections. Relevant order-a,, corrections to the expectation value of 
the four-fermion operator are drawn in Figs. 1 whereas Figs. 2 show the corrections to 
the formfactor. The gluon exchanges involving only light quarks merely renormalize 
the current in question, and we do not consider them. The corrections dressing the 
heavy-quark part vanish due to conservation of the 6-quark current (we consider 
gluon momenta much smaller than rrib)- 
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Figure 2: Diagrams contributing to the hght quark formfactor of a heavy hadron. 



In the nonrelativistic approximation for the light quark the "crossed" diagrams 
are suppressed, and the remaining diagrams Figs, la and b have obvious counter- 
parts in the corresponding diagrams in Fig. 2. Going beyond a simple potential 
approximation (e.g., at k'^ ^ ^q)y however brings in diagrams Figs. 1 c,d as well. 
In fact, one should keep in mind that in Eq. (pH) the integration of the formfactor 
is performed only over the spacelike components of q. This fixes the spacelike sep- 
aration of the bb and qTq currents to be zero, however per se does not specify the 
timelike separation of the vertices which is actually determined by the heavy quark 
propagators. In reality, a single diagram Fig. 2a corresponds to the sum of diagrams 
a and c in Fig. 1, and likewise with diagrams b. In the coordinate representation, 
the heavy quark propagator in Fig. la is ??(— Xq) and in Fig. Ic it is i^ixo) thus 
yielding unity in the sum (unity means absence of any propagation, as in Fig. 2). 

Let us illustrate it in the usual momentum representation. Denoting the gluon 
momentum in Fig. 1 by k, we keep the spacelike components of k fixed and consider 
the integral over uo = k^. The diagrams a and c are given, respectively, by 

Aik.uo] and A(k,uj] , (24) 

— — It ^ UJ — It ^ ' 

where A{k^u?j generically denotes the 'light' part of the diagram (including the 
gluon propagator). Since 

1 1 

\ = 2'Ki6{uj) , 

—oj — le uj — le 

the integration |^ of the sum of Figs, la and Ic amounts merely to setting fco = 
in the rest of the diagram. (This is a special case of the more general relations given 
in Appendix 1.) Then it exactly coincides with Fig. 2a if k is identified with the 
gluon momentum / in the latter. Although the gluon momentum transfer / is not 
generally equal to q but can differ by a primordial momentum in the bound state, 
integration over all g'is equivalent to integration over d'^/. Similarly, the sum of the 
diagrams in Figs, lb and Id yields the integral of Fig. 2b over q. 
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It is clear that this proof is generahzed for an arbitrary number of gluon ex- 
changes between the 'hght' and 'heavy' parts of the diagrams, or the case of the 
axial 6-quark current (see Appendix 1). It is imperative, however, that the b quark 
current is color-singlet. 



3 Applications 

We now turn to some applications of the relations (plj), (p2D. 

3.1 Perturbative renormalization of the color-straight oper- 
ators 

In general, the composite heavy-quark operators depend on the renormalization 
point fi which is assumed to satisfy the 'hybrid' hierarchy condition Aqcd <^ <C 
rrib- The most interesting is the logarithmic renormalization. This 'hybrid' renor- 
malization was first considered in ^, |^ where the one-loop hybrid anomalous 
dimensions were calculated for the quark bilinears and four-fermion operators. 

In the expressions of the matrix elements of the color-straight operators h{ak)hqTq 
via the integral of the transition matrix element of the light quarks current, the 
normalization-point dependence can appear in two ways: first, as a /x-dependence of 
the light-quark current itself. This is a usual, 'ultraviolet' renormalization since /i is 
an ultraviolet cutoff in respect to the light degrees of freedom. The second way the 
dependence on the UV cutoff can enter is via the divergence of the integral over the 
momentum of the final state. Indeed, in the effective theory with the cutoff /i the 
perturbative states with momenta above are absent, while the formfactors with 
|g I ^ /U coincide with those in full QCD. Therefore, if in full QCD the integral of the 



amplitude in Eqs. (pl| , |22| ) does not converge at g* ~ Aqcd but has a log behavior in 
the hybrid domain, this leads to the logarithmic dependence of the matrix element 
on 11. 

In practice we are interested in vector or axial currents of light quarks. They are 
conserved and their anomalous dimensions vanish (for the flavor- singlet axial-vector 
current there is an anomalous dimension in higher orders in related to the axial 
triangle anomaly). Therefore we will phrase our discussion neglecting this type of 
renormalization. 

The asymptotics of the actual light quark current formfactors of the heavy flavor 
hadrons is given by the perturbative diagrams where hard gluons transfer the high 
momentum from the light quark to the heavy one. The tree-level order-a^ diagrams 
are shown in Figs. 2a,b. By virtue of the relations Eqs. (0), ( P^ they determine 
one-loop renormalization of the four-fermion operators. It is easy to see that these 
diagrams yield amplitudes fading out at least as 1/g^ (the odd powers of q do not 
contribute to the integral). [| In principle, depending on the particular form of F, the 

^Similar quark counting rules in heavy mesons for \q \ <C vfib have been appUed, e.g., in ]lO| . 
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asymptotics of these diagrams may have the term - it is given by 

Ajccy 1 ~ — 

^^T^WWhq{^oiT-Ttto)t^q{^)\H,m ■ (25) 

{1)1% baf^b for the spin-flip transitions). The matrix element may not vanish 
for beauty hadrons with nonvanishing spin of hght degrees of freedom (let us recall 
that qo = 0). However, in this matrix element both hadrons are at rest, therefore 
any such 1/|^|^ term vanishes upon integrating over the direction of q. The fact 
of vanishing of the leading-order hybrid anomalous dimension for the operators of 
the form (bb) (^7^(75)5') was noted in [§] already in the mid 80's as a result of 
simple calculations of the one-loop diagrams. Our relation gives it an alternative 
interpretation. 




Figure 3: Order-a^ diagrams determining the asymptotics of the light quark form- 
factor of a heavy hadron. Similar diagrams with the twisted gluon lines are not 
shown. 

A closer look reveals, however, that the cancellation of the leading 1/|^|^ asymp- 
totics does not hold already at the one-loop level. The asymptotics has actually the 
form ~ ctK*?*)/!^!^ which emerges from the diagrams shown in Figs. 3 (other dia- 
grams decrease faster in the Feynman gauge). This leads to a nonzero anomalous 
dimension of the color-straight operators at order and their mixing with color- 
octet operators. In particular, the evaluation of the one-loop amplitudes leads to 

{Hb\bbqrq\Hk) ( . _ ±\ {Hb\bt% qTty\H,) ' 

(26) 

Eq. ( pH ) then yields for the UV part of the four-fermion operator 



■ In — - + finite piece , 
(27) 

and likewise for the spin-triplet operators. 



Ar{q) 



2 s 



4 



— j66gTg + NJl-—jbt%qTt'^q 
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A direct calculation of the two-loop anomalous dimensions confirms this. The 
computational details are described in Appendix 2. Here we only quote the result. 
Let us denote | 

Then 

7.1 = ''K (l - 7^) (l - 7^) (H) = 3iV. J + 0(ai) . (29) 

We note that since 711, 712 and 721 vanish to order as, these two- loop anomalous 
dimensions do not depend on the renormalization scheme. For 722 the second-order 
terms depend on the scheme and we do not consider them. 

Additional terms are present for the flavor- singlet operators: for the vector cur- 
rent only 722 is modified, 722 — > 722 ~ |^/f^- If the operator has the flavor- singlet 
axial current then only the diagonal anomalous dimension for the color-straight 

operator changes, 711 711 - 6nf {^N^ - (f^) . 

It is interesting that, although 711, 712 and 721 already appear in the second loop, 
they are universal. In particular, they are the same for both timelike and spacelike 
components of the light quark currents. A priori this does not need to hold. We 
expect that this universality will be violated in the next order in as- 

The two-loop anomalous dimensions are enhanced, they contain a large factor 
7r^. Neglecting them introduces a numerical uncertainty in the running of operators. 
We can estimate it by simply setting In ^ to unity. The corresponding corrections 
at as = 1 constitute about 15 to 30%. This provides additional justification for the 
standard choice of as(/i) = 1 as the low (hadronic) normalization scale. 

We point out that the naive estimate of the power of the asymptotics of 
the light current formfactors existing in the literature is not correct: the actual 



fall off is only 1/|^|'^ as shown in Eq. (|26|) , which, however, is generated only by the 
exchange of two gluons with momenta ~ q. (The modification for the spin-triplet 
operators is obvious). This asymptotics can be easily RG improved using relations 
Eq. (^), (^. To the NLO it amounts to adding the factor [as(g)/as(/u)]"^^'^^^° 
(or [as{q) / as{fi)]^~^^''^^^^'^^^^'^^° for the fiavor-singlet vector current) in front of the 
color-octet operator, with /i the normalization point of the operators. Since the gap 
between the typical hadronic mass //had and is not too large in the logarithmic 
scale, this observation has, probably, rather theoretical than practical significance. 

■^The anomalous dimensions of the operators are often defined with the opposite sign. We prefer 
to use this convention where the meaning of the anomalous and canonical dimensions are the same. 
That is, the scaling properties of the operators are given by the sum (rather than difference) of 
their canonical and anomalous dimensions. 
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It turns out that in the relativistic system the Coulomb interaction is still strong 
enough to make the wavefunction density not approaching a literal constant at zero 
separation but having the small logarithmic dependence on distance which appears 
only at the level of loop corrections, at order ~ a^. 

Based on the application to the lifetimes of heavy hadrons (first of all, in B 
mesons) and routine application of factorization, a standard choice for the basis for 



four-fermion operators ascending to the original papers on the subject was 11, 12 



Osingi = (&r(i)g) {qT^%) , Oo,t = {bf^T^'k) {q^T^^) (30) 

(t"' = A" are the usual Gell-Mann color matrices), that is the s-channel color- 
singlet and color-octet operators. It appears, however, that a better choice is to 
classify the operators according to the color structure in the t-channel: 

o = {kr^'k') (gfcr(2)6') . 5uSkj , t = {kr^'k') {qkT^'^b') ■ tTkj ■ (3i) 

In the large- Ai'c limit these two bases coincide (up to permutation): 

0.,.„ = 2T + i-0. 0„„ = 1 (l - ij) O - i-r (32) 

O = 20oct + -^Osingl , T = - ^1 - — ^ Osingl - ^^oct • (33) 

The t-channel octet operators T also diagonalize the one-loop anomalous dimension 
matrix; its value depends on the type of the current, fiavor-singlet or octet [Q. 

We parametrize these generic expectation values encountered in actual weak 
decays as 

^ {B\Ov\B) = uv ^{B\Tv\B) = Ty (34) 



1 



(B\Oa\B) = UA ^{B\Ta\B) = TA . (35) 



The parameters u, r have dimension and are constants in the heavy quark limit. 
They can be valence or non-valence; the flavor of the light quark in the operator 
will be indicated as a superscript. 
For Af,-baryons we denote 

■{Ab\Ov\Ab) = A -J—{A,\Ty\A,) = -h' ■ (36) 



2Ma/ " 2Ma/ 3 

in the valence approximation A' = A. The values of A for u and d quarks are equal, 
likewise for A'. These valence expectation values will be normally used without 
flavor index. 
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A remark is appropriate to conclude the discussion of the perturbative renor- 
mahzation. Strictly speaking, the flavor-singlet operators can be renormalized in 
somewhat different ways depending on the prescription to treat the tadpole-type 
closed loops. The free quark loop by dimensional counting scales with the UV cutoff 
H like /i^, and describes a possible power mixing with the D = 3 "unit" heavy-quark 
operator bb already at order a^. Although for practically relevant operators such a 
"bare" mixing vanishes for the usual way to regulate the light quark loop, one can 
raise the question where this freedom is reflected in relations (21,22) for a generic F. 
The resolution is rather straightforward: the flavor- singlet current qTq also requires 
regularization of the closed fermion loop and, a priori admits mixing with the unit 
operator (the tadpole graph). This operator does not lead to any physical transi- 
tion at g 7^ but to the forward amplitude with g = 0. A formally defined current 
qTq may thus lead to an additional term proportional to S^{q) in the transition 
amplitude A{q), which would reproduce the tadpole term in the expectation value. 

Similarly, strictly speaking one could have chosen an arbitrary convention for the 
phases of the states |iff,(g )) with different momenta q. This would redefine the phase 
of the transition amplitude A{q). In our relations such a freedom was eliminated by 
adopting Eqs. (|TT],|T2],|Tip which ensures, for example, the proper analytic properties 
of the transition amplitudes. 

In the purely perturbative calculations one can, in principle, consider not only 
the actual physical amplitudes, but also similar transition amplitude induced by the 
light quark currents carrying color. Applying to them relations similar to Eq. ( |2TD 
and (^) one would need to consider the color-nonsinglet quark in or out states. 
This case requires certain care since such amplitudes may have additional (gauge- 
dependent) infrared singularities. 

It is worth reiterating that in our analysis it is assumed that all heavy quark 
operators are renormalized at a scale well below rrih, which implies a nontrivial - 
even if finite - renormalization when passing from the full QCD fields. In particular, 
the vector bb and axial bab currents both do not renormalize in this domain; however, 
they run differently when evolved down from the scale ~ rrib. While 6706 — > bb 
is not renormalized, the short-distance renormalization of Irfklhb — ^ bakb slightly 
suppresses it: 

(the second-order correction has been also calculated [|I^]). This is not, however, 
the only short-distance contribution differentiating the renormalization which, in 
general, depends on the exact form of the operators. We will not further dwell on 
these corrections in our numerical analysis. 
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3.2 Estimates of the color-straight expectation values in B 
mesons 

Relations (|l],|2D open a possibility for an alternative evaluation of the expectation 



values of the color-straight operators. It requires knowledge of the light-quark- 
current formfactors of heavy hadrons. The direct experimental information about 
them is scarce. Therefore, we have to assume a reasonable model. Our general 
strategy for all expectation values of interest is the same: decompose the transition 
amplitude into the invariant formfactors, and adopt a model for the formfactors 
satisfying known constraints. 

For the family of A;, baryons, the number of possible amplitudes is limited due 
to the fact that the light degrees of freedom are spinless - one can construct only 
scalar, vector and tensor currents while pseudoscalar and axial amplitudes vanish. 
There are no axial analogues of the expectation values in Eq. (|36|). For mesons all 
amplitudes are possible. Our main attention will be devoted to the vector and axial 
currents, due to the chiral invariance of phenomenologically relevant four-fermion 
operators. We do not consider the tensor current, and only briefly comment on the 
scalar one. 

There is only one formf actor for the vector current (for each flavor content) for 
both B and Af, describing the only nonvanishing timelike component: 

{B{q)\J,\BiO)) = -v.Fsiq') 

(Afe(g)|JjA,(0)) = v^Fj,^{q^)u{v,s')u{v,s) . (37) 

One important constraint on the formfactors is their value at = 0. The values of 
^B,At (0) s-re fixed by the corresponding charge of the hadron: it is 1 for the current 
of a valence quark, and zero for a 'sea' light flavor .0 For the amplitude in Eqs. 
the integration over q yields 



{2tt 

The valence formfactors are expected to decrease for non-zero q^. For the isovec- 
tor formfactor the slope at = (related to the corresponding charge radius) can 
be estimated in terms of experimentally observable quantities by an analogue of the 



Cabibbo-Radicati sum rule for heavy hadrons ||I4 



^ -^EP^+D^^^^t^ (39) 



where the sum runs over excitations of B with spin J and Q is the light quark 
charge in units of e (a similar relation holds for baryons as well), and the nonresonant 
contributions are neglected. This slope is not yet known experimentally well enough. 



^We adopt the convention where B mesons have the quark content bq. 
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Relations following from Eq. ( ^Sf ) can be used in the lattice simulations to eval- 
uate the expectation values, by measuring the transition formfactors in a few kine- 
matic points and interpolating between them. This type of lattice measurements 
can be simpler than for the heavy-quark current transitions, since the heavy quarks 
remain at rest and the momenta involved in the process do not scale with mq. This 
makes the static approximation rather straightforward. 

If we represent the formfactor as a sum over singularities in the t-channel 

n^) - E ^ . (40) 

the integral Eq. (^) takes the form 

i^,nf)--i^c.M!-, (41) 
we have an additional constraint 

E CnM^ = (42) 

n 

following from the fact that the transition amplitudes decrease faster than 

It is natural to consider the simplest model of saturation containing only two 
lowest-lying 1~ states with appropriate isospin quantum numbers. For example, 
for / = 1 we use p(770) and p(1450). With fixed normalization at = and 
the constraint ( ^2]) this model predicts the value of the integral in terms of the 
two masses. It is worth noting that such a model would obviously lead to equal 
expectation values of the operators in B and A;,. Imposing an additional constraint 
from the slope of the formfactors would allow one to fix all residues in a three-pole 
model as well, which can be hoped to yield a more accurate estimate. 

A word of reservation is in order at this point. Such a saturation of the nucleon 
formfactors by two lowest t-channel resonances is known to provide a good approx- 
imation for moderate where the experimental formfactors are described by the 
double-pole expressions. There is no general theoretical justification for such a coin- 
cidence, and more resonances are expected to play a role for larger q^. In particular, 
at — > 1 GeV^ the formfactor can decrease faster. Due to the phase space factor 
the role of the domain of large is enhanced. The contribution of higher states, 
while affecting a little the formfactors near g^ = 0, still can significantly change the 
integral (^T]). We will return to this point later. 

It is clear that since the asymptotics of the amplitudes has an odd power of 
l/|g|, their representation by a finite number of the t-channel resonances is not 
possible. The true spectrum of the t-channel states must extend to arbitrary high 
masses. It applies even if there were no typical for QCD log-like dependence of the 
asymptotics l/|g |^. It does not affect our estimates since we evaluate the operators 
in the effective theory where the high-momentum component of the hadrons is peeled 
off. 
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Addressing the color-straight operators containing the axial light-quark current 
(which does not vanish in B mesons) we note that its matrix elements are generally 
described by two formfactors, just as in the well-known case of spin-i fermions. 
These are analogues of the axial-charge and weak magnetism terms. Spontaneous 
breaking of the chiral symmetry modifies the value of the axial-charge formfactor 
from its symmetric limit of 1 at = 0. Nevertheless, for the isovector current, 
its conservation d^J^^{x) = in the chiral limit leads to a relation between the 
formfactors, so that only one, the axial-charge formfactor is independent, as in 
the case of the vector current. At = this relation equates the axial-charge 
formfactor to the B*Bn coupling g (the heavy-quark analogue of the Goldberger- 
Treiman relation). Given the value of g, therefore, one can evaluate the expectation 
value exactly as outlined for the vector currents. 

For the isosinglet axial current, one has to take into consideration the anomalous 
term, the topological charge density Q: 

d.jjSix) = 2tqmg^,q{x) + n^Qix) , Q{x) = ^ Tr G,^G°^(a;) (43) 

with ihg the light quark mass matrix. The matrix elements of Q over the B meson 
states are not known, and the above relation appears to be less constraining. In the 
large- A'^c limit the difference between singlet and nonsinglet formfactors is expected 
to disappear; however, the practical validity of this approximation for the anomalous 
term is questionable. These problems are addressed in the next section. 

Let us briefly mention the case of the scalar current. Although the corresponding 
formfactor is not fixed at = 0, its value for the valence quarks can be obtained 
from the SU (3) mass splittings: 

{B+ uu{0) B+) ^ ^ ^ 0.7, A, n« A5 ^ ^ 1.4, 

2Mb ms 2Ma^ 

(44) 

This estimate is obtained with the help of the Zweig rule in a similar way as 
done in to extract the value of the nucleon cx-term from the SU{3) splittings 
in the baryon octet. We neglected here the light quark masses mu,d and took 
ms(lGeV) ~ 130 MeV. The mass of the baryon Sf, has not been yet measured; 
the above estimate used the prediction M^^ = 5805.7 ± 8.1 MeV |]T6|. The normal- 
ization point dependence of rris in these relations reproduces the dependence of the 
scalar current. 

In what follows we will apply the described strategy to evaluation of a few ex- 
pectation values of operators with the vector and axial light quark currents. First, 
however, we discuss qualitatively the saturation of the integral of the formfactors in 
the adopted models. 

In the case of the (valence) vector current we have F{0) = 1. It is natural to 
think also that |F(g^)| < 1 at spacelike q. Let us further assume that F{q^) is small 
enough above a certain scale /i, so that we can neglect it there: 

F(g2) ~ at - g2 > /i^ . (45) 
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Then we get an upper bound 



2M 



B 



-{Hb\b-f^bq-ff,q\Hb) 



< 



0.017 GeV^ 



(46) 



for /i = 1 GeV. This bound is of the type discussed in 0; the numerical coefficient 
coincides with the one given there. 

Are assumptions hke Eq. (|^) reasonable? In the effective theory with the nor- 
malization point /i the momenta of fields exceeding /i are absent, whether or not 
the full theory yields a logarithmic 'tail' at large momenta. For example, it is not 
possible to exchange a gluon with momentum |g | > /i in such a theory. The exact 
shape of the formfactor would depend on the concrete realization of the effective 
theory. The amplitude may not vanish exactly due to multiple gluon exchanges 
with |g I < /X, however would then decrease exponentially. 

On the other hand, a step-like formfactor saturating the bound (^6]) is clearly 
unrealistic. Therefore, we can assume instead that 



< e 



(47) 



which results in 



2M 



B 



-{Hb\b-ff,bq-ff,q\Hb) 



< 



3/2 



0.022 GeV^ 



(4J 



with the same value for fi as in 



As a matter of fact, an exponential ansatz 
for the formfactor with fi^ adjusted to reproduce the 'charge' radius, is a 
reasonable model for the possible behavior of the valence formfactor of purely soft 
degrees of freedom. In particular, if we adopt the slope at = following from 
the pole model keeping only the two lowest states with masses Mi and M2 (this 
is a good approximation in known cases) and use the exponential formfactor, the 
resulting integrals appear noticeably smaller than in the two-pole ansatz itself: 



1 



2M, 



B 



■{Hb\h,bqj^q\Hb) = F(0)- 



1 



SVT' 



3/2 



' MfMj 
Ml + M| 



3/2 



(49) 



MfMl 



The above discussed bounds rely on the assumption that |-F(g^)| < 1 . It always 
holds in nonrelativistic QM, however we do not know a general rigorous proof in 
QCD. For the isovector current one can employ the equal-time commutation relation 
(J^ can be both the vector = q'Jfj.^r'^q or the axial current = g7^75|T'*g) 



J+{x),J^{y)] = 6'{x-y)2V,\x) 



(50) 
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to represent 1 — |-F(g^) p at < as a difference of two sums of tlie distinct transition 
probabilities: 

= 1 - (E \Fn? - E \F^n?] ■ (51) 
\ n m / 

Here l-^mP scliematically denote tlie transition probabilities in, say, B~ meson 

induced by the currents U'-y^d and d'-you with the momentum transfer g, respectively 
(and similarly for A;,). In the second sum only the states with / = | contribute. 
Since there are no valence d quarks in 5", in the large- iVc limit the last term with 
the wrong sign would vanish. Also, in this limit the isoscalar meson formfactor 
is expected to coincide with the isovector one. Therefore, the large-A^^c arguments 
allow to establish such a QM bound for all formfactors of interest. 

There is no natural normalization for the axial (pseudoscalar) formfactors at 
small momentum. Moreover, the amplitudes generally have an enhancement due 
to the pion pole. However, the domain g*^ ~ yields a very small contribution 
to the integral (see, e.g., Eq. (^61) ). The significant contribution can originate only 
from momenta > 1 GeV where one expects the effects of chiral symmetry breaking 
to become insignificant. As noted above, the equal-time commutation relation (|50[) 
can be used to derive a sum rule of the type (^TJ) also for the matrix elements of the 
axial isovector current. Its explicit form is similar to (|51|) and reads 

\G.{f)? = 1 - (e l^nP - E \G7A (52) 

with Gi(g^) defined below in (|68| ) and IC^P, jCmP are the analogues of the F„ 
amplitudes for transitions induced by the axial current acting on a S meson. At = 
this sum rule is just the familiar Adler-Weisberger sum rule and the amplitudes Gn 
are related to pion couplings between the ground and excited states. The explicit 
form of these sum rules for heavy mesons and baryons can be found in |]14|, p!7| . 
Therefore, we expect the type of bounds (4^,4^) to hold also for the axial current 
expectation values as well. 



4 Numerical estimates 

In this section we estimate the expectation values of the color-straight four-fermion 
operators relevant for the lifetimes of beauty hadrons. The light quark fields are left- 
handed; the Penguin diagrams bring in the right-handed fields as well. Nevertheless, 
the chiral structure of the currents admits only the vector or axial light quark cur- 
rents. Since the coefficient functions can include the momentum of the decaying b 
hadron (its velocity), the timelike and spacelike components enter, in general, with 
different weights; the three-dimensional rotation invariance is still preserved. Fi- 
nally, since the forward matrix elements are considered, only the parity-conserving 
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(three-dimensional) scalar expectation values survive. Therefore, we need to con- 
sider the operators 



Ov = {bb){q^oq) = ip-i^h){q-i^q) , Oa = -{bab){q^^5q) = (&7M75&)(g7M75g) • 

(53) 

As was mentioned, the expectation value of the operator Oa in vanishes. 
4.1 Vector current 

We first consider the case of B mesons. Assuming only isospin symmetry, we define 
the isovector and isoscalar four-quark matrix elements by 

-^{B,\b^^bqT'^^,q\B,) = VsT^^ 



2Mb 
1 



2Mb 



mh,b 11,l\Bj) = V,5,, (54) 

q=u,d 



The indices i,j label the respective state in the isospin doublet i = {d, —u). Accord- 
ingly, we introduce the isospin-triplet and singlet vector formfactors 

^ (5.(g)|gr"7Mg|5,(0)) = v,:F,{q^)Tt^ 



2Mb 

1 



^ q=u,d 

with the normalization conditions J^i(O) = ^^3(0) = 1. Using the two-pole ansatz 
saturated by p(770) and p(1450) for the nonsinglet current JF3 , we get from Eqs. (^ 

H 

1 M'^ /lf2 

^ _J_ M.ivi, ^ 0.045 GeV^ . (56) 
47rMi + M2 ^ ' 

It is natural to saturate the / = formfactor JF^ by the states u{782) and u;(1420). 
It then leads to almost the same numerical estimate for Vi as for V3. The reason 
is obviously an almost exact degeneracy of the vector states in the isovector and 
isosinglet channels. Although it perfectly fits the large-A^^c picture, we cannot be 
sure what is the actual accuracy of such a conclusion. Nevertheless, in view of 
such a suppression of the difference, we will take the two different combinations 
of the expectation values which actually parametrize the valence and non-valence 
contributions: 

"^^y^ = ^{B-\h,buj,u\B-) ^ -0.044GeV^ (57) 



2M, 



l^{B-\h,bd^,d\B-) ^ O(l0-^GeV=^) . 
vib ^ ^ 
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The last number, clearly, is at best an order of magnitude estimate. 

We can try to estimate the violation of the SU{3) flavor symmetry considering 
the expectation value of fry^bs'yfj^s in Bg mesons. For this we saturate the formfactor 
with the vector ss states 0(1020) and 0(1680), which corresponds to the "ideal" 
mixing in the uj — (j) system [|1^ . In this case we would get 

^ (5, 167^6 S7^s|5,) ~ -0.085 GeV^ (58) 



2M, 



i.e. almost twice larger than the first estimate (p7|). 

A closer look reveals, however, that the above expectation values are saturated 
at rather high momenta. Half of the 'valence' value comes from |g | > 1.5 GeV, and 
from even higher momenta in Bg. For this reason these estimates exceed the bounds 
(1^), (^8|) discussed in the previous section, for a reasonable scale /i ~ 1 GeV. For 
example, adopting the exponential ansatz for the formfactor, we would get 

-^{B-\h^.hu^^u\B-) ~ -0.007GeV3 

^{Bs\h,,bsj^s\Bs) ~ -0.015GeV^ (59) 

A somewhat unexpected result of these simple estimates is the apparently large 
amount of SU{3) breaking in Eqs. (0,|^). While it is not clear to what extent this 
is an artifact of our use of the simple two-pole ansatz for the formfactors over a wide 
domain of q^, it is worth noting that a simple mechanism exists which could account 
for it. It is well-known that the isovector charge radius of a hadron diverges in the 
chiral limit |T9. This indicates that the contribution of the low-momentum region 



in the integral over the formfactor (|38D is more suppressed in nonstrange B mesons 
compared to the case. Since the two-pole model does not capture the origin of 
this phenomenon (the contribution of the two-body vrvr intermediate state in the 
t-channel), it is conceivable that the magnitude of SU{3) violation in the matrix 
elements does exceed a few percent. 

It is interesting to compare the above estimates with the evaluation based on 
vacuum factorization. Both types of estimates have the same O (iV°) scaling in A^^- 
However, the vacuum contribution for the color-straight operators is O (A"°), similar 
to other meson states. This is in contrast to the case of the t-channel octet operators 
where the vacuum state is A'c-enhanced. Therefore, the factorization estimate is not 
expected to give an accurate result. For the valence expectation value one has 

1 ,^ , _ 1 



{B-\h^bu^,u\B-)i,,,,, = -^fl{^)MB (60) 



2Mb 

(the non- valence value vanishes). Here Jb denotes the annihilation constant of B for 
the h'^al^u current normalized at a low point ^ where factorization must be applied 
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0, H, in contrast to the physical fs defined for the current normahzed at /i ^ mj,: 



fBifi) ~ Jb —, r • 61 

[_as[mb)_ 

The physical value of lies, probably, around 160 MeV. However, to the leading 
order in l/m^ we work in, it is more consistent to use the asymptotic value which 
differs from the physical one by 1/mf, and non-logarithmic perturbative corrections. 
These decrease the physical value of fs by about 20% Therefore, we adopt 

fs = 160 MeV for q;s(/x) = 1, yielding 

^ (E- 167^6 n7^n|5-)factor ^ -0.011 GeV^ , (62) 



2M, 



which is significantly lower than Eq. 

The fact that the corrections to factorization can be significant, is expected. 
Unfortunately, there are good reasons to question the accuracy of the alternative 
estimate (0) either, and a too large SU (3) breaking is another indication. We think 



that it is justified to consider the estimate (|57|) for the valence expectation value 
rather as an upper bound, while the number obtained in the exponential ansatz a 
reasonable lower bound. A conservative estimate then is 

^ {B-\h^bu-f^u\B-) = -(0.025 ± 0.015) GeV^ 



2M 



B 



-L-{B-\h,bd^,d\B-) ^ O (5 ■ 10-^ GeV^) . (63) 

Similar estimates can be adopted for strange quarks in Bs- 

Next we turn to baryons. Under the light fiavor SU (3) group the A;, and states 
transform as an antitriplet Tj = (H^, — H^, A^). In the limit of SU{3) symmetry there 
are only two independent formfactors, which can be defined as 

■{T,{q)\qX^^,q\T,{0)) = v, J^^{q') \%u{v, s')u{v, s) 



2Ma, 
1 



2M ^^^^^"^l ^ l^f^ll^m = v,Tt{q')6.,,u{v,s')u{v,s). (64) 

q=u,d,s 

The normalization at = is ^^(0) = 2, J^^(O) = —1. Using a similar model 
for the formfactors as in the meson case we get the same expectation values (up 
to the sign) for the valence matrix elements, and strongly suppressed non-valence 
contributions. For example, the two-pole model yields the following value for the Af, 
matrix elements 

1 , . i7 , _ I. \ 1 / . i7 , 7 ,1 . \ I 0.007 GeV'^ (exponential) 
^(A,|67.&«7.«|A,) = ^(A,|67M&rf7.ci|A,) ^ | o.045GeV3 (two-pole) 

(65) 
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For the same reasons as before it is natural to consider the two-pole value as an 
upper bound. The expectation values of the non-strange operators in the states 
emerge the same as in (^), whereas (Sfe|67^6 s7^s|Sfe) again hterally appears twice 
larger than in the SU{3) limit. As discussed in the B meson case, such a large 
symmetry violation can be suspected to be, at least partially, an artifact of the 
two-pole model. 

It is worth noting that in the case of heavy baryons the light quark formfactors 
have anomalous Landau thresholds associated with the NNB triangle diagrams. It 
is well known that it is such singularities that determine the low-momentum behavior 
of the formfactors and, in particular, the large charge radius of weakly-bound states 
like deuteron For the formfactor the anomalous singularity starts at 

In this system, however, the corresponding mass still lies higher than the states we 
use to saturate the formfactors. Moreover, there is no reason to expect the residues 
to be significant (for example, they are 1/A^c suppressed). Therefore, we believe that 
these singularities do not play a role in the expectation values we study. In any case, 
a refined estimate will be possible with a better knowledge of the formfactors, say 
using determination of its slope based on the application of the Cabibbo-Radicati 
sum rule to the radiative decays of excited baryons. 



4.2 Operators with axial current 

The expectation value of the operators &7/x75&^7/x75'? vanishes in the A;, baryon 
family, and we consider it only for B mesons employing the relation Eq. (|2^). In 
this case 

(Sbe)|S(g)) = |S*(g,e)), = ^ d'xb^^.bix) . (67) 

Since the light degrees of freedom carry spin |, the axial current is parametrized by 
two formfactors; the third possible structure has wrong T parity and vanishes. This 
is an exact analogue of the absence of the second-class currents in /^-decays of light 
baryons. Thus one has 

^{B:{q,e)\ E g7,759(0)|5,(0)) = {e;G'f\q') - {e*q)q,G^^\q')} 6., , 

^ q=u,d,s 

-L_{B*{q,e)\qX^^,^,qmB,{0)) = {e;G,{q') - {e*q)q,Go{q')} X^, . (68) 

Absence of the structure {€*q)v^ is easy to show explicitly (note that in any case the 
timelike component of the axial current does not enter the four-fermion operators). 
Using Eq. (^) and the fact that Sb commutes with all light-quark field operators, 
we get an equality 

(i?*(g,6)|J,5(0)|i?(0))* = (i?(0)|J,5(0)|i?*(g,6)) = (i?*(0, e^)| J,5(0)|i?(g)) . 
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Inserting here the formfactor decompositions (|68|) for these matrix elements and 
taking into account the fact that Gi are real from T invariance, one finds that the 
structure (eg)f^ appears with opposite signs on the two sides of the equality. Hence 
its coefficient must vanish. 
We thus get 

1 — — 1 

■{B-\h^^,b{u^^^,u-dj^j,d)\B-) = -— / dtVi{3G,i-t)+tGo{-t)) 



2Mb ' '^'-^ ^ "^'^ /I / 

(69) 

and a similar expression for the singlet matrix elements. 

In the chiral limit, which will be assumed in what follows, the isovector formfactor 
Gi at = is related to the BB*tt coupling: 

Gi(0) = g . (70) 

The nonrelativistic quark model predicts g = —0.75. However, the QCD sum rules 
estimates yield lower values |^ |2^, The recent analyses predict g = —0.3 
|, 1^ which is consistent with the existing experimental bounds g^ = 0.09 — 0.5 
. Moreover, the equation of motion d^J^^ = leads to q^G^^q^) = Gi{q^) at all 
g^, therefore for the nonsinglet expectation value Eq. (|69D takes the form 



1 — — 1 

{B-\h^-1,h{u^^-i^u-d^^^,d)\B-) = -— / dtv^Gi(-t). (71) 

The only nonvanishing contribution to Gq from the pseudoscalar states in the 
isovector channel comes from the massless pion. The J^^ = states contribute 
to both Gi and Go- 



n 



Goit) = + E (72) 

ml-t ^ Ml-t 

with the condition 9n = 9 replacing the zero-transfer normalization of the vector 
formfactor. A faster than fall-off of the transition amplitude requires addition- 
ally 

E 9nMl = , (73) 



which is analogous to the second constraint in Eq. (^21) for the vector current. 

In the numerical estimates for the isotriplet current we will consider both a 
two-pole ansatz for Gi(g^) and the exponential ansatz 

G,(g^)=^e-V-\ Go(g^)=^^ with l^' = ■ (74) 

Such a choice of fi ensures that the two ansatze have the same behavior at small g^. 



23 



In the I,J^'-' = 1,1^^ channel only the lowest-lying state ai(1260) has been 
observed. For the numerical estimates we will need also the mass of its first radial 



excitation a[. This has been extracted in [27| from an analysis of the Weinberg sum 



rules. The value obtained in for the mass of the a'l resonance is 1869 MeV, 
which is what we will use in our estimates. 
First, with the two-pole ansatz we obtain 

1 - - a M'^M'^ 

-{B-\h,^,b{u^^^,u - d^,^,d)\B-) = -A^^ ^ 0.084GeV3. (75) 



2Mb' 27rMi + M2 



As explained above, we have adopted in this estimate the value g = —0.3 
For the exponential formfactor one obtains a smaller value 

(76) 

In a completely analogous way one can estimate the matrix element of the / = 
octet axial current with the flavor content of rj. With the mass of the state /i(1285) 
close to mass of ai and assuming a similar degeneracy for the second excitation we 
do not get appreciable SU (3) violation and, therefore, obtain for 



^(5 |(fe7M75&)(M7M75M + rf7M75C?-2s7^75s)|5 ) 



2M, 

the same value as in Eqs. ([76|) and (ffSl). 

In the case of the singlet axial current we need to account for the presence of the 
anomalous term in its divergence, 

d,jf^{x) = UfQix), Q{x) = ^GG{x). 

For simplicity, we will assume the exact SU{3) chiral limit. The value of the isos- 
inglet axial formfactor at small is given by the matrix element of the anomalous 
divergence Q{x): 

^ {B*{e,q)\QmBm = -Gf\o){e*q) + 0{q') . (77) 



We used here the fact that there are no massless particles in the singlet channel and 
consequently C'^^ is finite. The contribution of the pseudoscalar states to the form- 
factor Gq'^ does not vanish and is determined by their coupling to Q{x), {n\Q{0)\0). 
Similarly, the Gf*^ and Gq'^ formfactors at arbitrary are not directly related to 



each other, but the difference i{Gi'^ —q'^G'q^) equals to the matrix element oinfQ{ 
Very little is known directly about these fiavor-singlet expectation values or B*Br]^'^ 
coupling. 

Nevertheless, for estimates one usually employs an approximation in which the 
matrix elements of Q{x) are saturated by the r]' pole. Moreover, the couplings of the 



x) 
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whole nonet of the pseudoscalar mesons vr, t], rj' are assumed SU (3)-symmetric. 
This assumption is incorporated in the simple cr-models proved to be successful 
in describing the properties of light hadrons. This model ^ naively has an 



U{3) X U{3) chiral symmetry; the U{1) problem is solved by adding the anomalous 
term with Q{x) and assuming the nonvanishing (in the quenched approximation, 
that is, in QCD without light flavors) value of the zero-momentum correlator of the 
topological charge densities Q{x) 

= J d'x{0\zT{Q{x)Q{0)}\0)n,=o (78) 

which leads, basically, to the nonzero anomalous mass of t]' meson m"^, = A^//^. 

Adopting such a model, we also have Gi'\o) ~ Gi(0). 

A possible justification for such a picture lies in the large-A^^c approximation. 
However, in this limit m^, oc 1/Nc and the anomalous U{1) symmetry effectively 
restores, which seems not be close to actual world where the anomalous mass of 
77' is numerically large and the octet-singlet mixing in pseudoscalars is small. It is 
probable that there exists a deeper dynamic reason explaining the practical validity 
of such approximation. 

The model with a single rj' state in the pseudoscalar channel which merely shifts 
the pole in the nonsinglet amplitudes from = to m^,, while describing reasonably 
well the low-g^ matrix elements of the topological charge density, leads to their too 
mild suppression at large g^. In reality they are expected to decrease very fast above 
a typical momentum scale of the nonperturbative vacuum configurations. In order to 
mimic this behavior, we have to employ at least two pseudoscalar states saturating 
the correlators of Q, and we take the state t]' with a mass of M^/ = 1295 MeV as 
the second pole. One expects an J^'^ = SU{3) singlet in this mass region, 
accompanying the observed octet of pseudoscalars containing 7r(1300), ?7(1295). In 
reality, the wide 'gluonium' states can give a significant contribution. Probably, an 
exponential ansatz is a better approximation here. 

In principle, the spectrum of the axial- vector singlet states has no direct relation 
to the anomaly and the U{1) problem. Hence we take for the corresponding masses 
the experimental values, namely /i(1285) and its first radial excitation f[, neglecting 
their mixing with the octet states. /i(1285) lies close to the isotriplet state ai(1260), 
indicating smallness of the annihilation effects. Therefore we will take for the mass 
of the first radial excitation Mj^ = Ma'^ ~ 1870 MeV in the numerical estimates 
below. In the two-pole model we have 

G(0)(g2) = ^(^^(O) 



{Ml-q^){M},-q 

W'-?^)K'(i295)-? 
The last equation replaces the second of Eqs. (|72| 
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As a result, the difference in the estimates compared to the isotriplet current hes 
basically in the anomalous term, and is not too significant. Numerically, we get for 
the two-pole ansatz 



~ -(0.29 + 0.05A)Gf\0) GeV^ = 0.1 GeV^* at Gf\o) = -0.3 . (80) 

In the numerical estimate we used the equality Gf'\o) = Gi{0) = g which holds in 
the large- iVc limit, as discussed above. 

We present also a calculation of the matrix element (^) employing the exponen- 
tial ansatz. This is constructed in the same way as for the axial charge formfactor. 
For the topological charge formfactor we use an exponential normalized at = by 
the same value Gf'\o) and vary the slope parameter Hq from fig = m^, = 0.92 GeV^ 
(corresponding to the pole dominance by r]' alone) to 0.59 GeV^ (corresponding to 
the two-pole model, see Eq. (^9])). This yields the following numerical estimate: 



1 .o-,i. . . j2 - - - -- - 

q=u,d,! 



/r 



Ml + MJ, 



~ -(0.053 + (0.010 to 0.020))G'!°^(0) GeV^ - (0.02 to 0.023) GeV^ (81) 

with the same value for Gi'\o) as before. One could try to estimate the effects 
of SU{3) breaking by accounting for the known shifts in masses and mixing. We 
think, however, that such models are too crude to capture correctly details of SU{3) 
violation, and we do not attempt it here. 

Combining the above results for the octet and singlet expectation values, we get 
the following estimates for the valence and non-valence axial expectation values: 

^ 'B-\}ri^'-l^hu'^^-i^u\B-) ~ 0.018 GeV^ (82) 



2M, 



2M 



B 



0.09 GeV^ 

(fi- 1 67^756^7^75^1 5") - 0.002 GeV^ (83) 

0.007 GeV^ 



where the upper (lower) value corresponds to the exponential (two-pole) formfactor 
model. 

We note that the effect of the axial anomaly can be numerically important, 
although it is formally of order l/iV^. In the approximations considered here, it 
makes a contribution of about 15% of the total singlet expectation values (|80D and 
(pl|), respectively, and it can dominate the non- valence matrix elements. 
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Finally, we quote also the factorization approximation estimate for the valence 
expectation value of the axial current. We obtain 

^(S-|67^756M7^75M|S-)factor = ^/|(/x)Mb ~ 0.034GeV^ . (84) 

Just as for the color-straight operators containing the vector current, we do not 
expect the factorization approximation to be accurate. However, it is interesting 
that the factorization value for the axial operators is less suppressed compared to 
the case of the vector current, and appears to be closer to the estimates given above. 



4.3 Estimates from the fourth sum rule 

One of the color-octet operators, the flavor- singlet vector four-fermion operator can 
be estimated in an alternative way. This operator Od 

Od = ih,t%)iq7,t^q) (85) 

q=u,d,s 

is the QCD generalization of the Darwin term in atomic physics ^, 

-J—{H,\27rasdn\H,) = --^{H,\On\H,) = - {pl) ^ ■ (86) 

On the other hand, it determines the third moment of the small velocity (SV) struc- 
ture function, the so-called fourth sum rule, and is related to quantities measurable 
in the semileptonic decays. For example, in B mesons this sum rule in the resonant 
approximation takes the form |R2[ 



IpI = ^V2ki/2|' + 2Eyn/2\' + ■■■ , (87) 

where Tj are so-called "oscillator strengths" which determine the small velocity tran- 
sition amplitudes into the excited p-wave states with spin of light degrees of freedom 
j. The excitation energies of these states with respect to the ground state s-wave 
mesons are denoted with Ej (for a recent discussion see review [^). 

It should be noted that the literal application of the fourth sum rule requires 
specific regularization scheme for the operators. In view of the tentative nature 
of our estimates we neglect these subtleties here. Some of them were discussed in 
IP and more recently in and |^. We only mention that the large negative 



logarithmic anomalous dimension of the Darwin operator 

(asOn),' ^ ( ^ j {a^Oo), , /3o = yiVe - -Uf (88) 

to a large extent offsets the apparent scale dependence in evaluation of (Od) via p% 
due to the //-dependence of (or, similarly, the scale dependence of the factorization 
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estimate of p|) in B mesons). According to the standard practice we use for our 
estimates the scale corresponding to as{fi) = 1. 

The recent discussion of the status of the sum rule evaluation of pf^ in B mesons 
can be found in the review Sect. 4. The corresponding value is in a reasonable 



agreement with the factorization estimates (see Appendix 3, Eq. ( |A3.7| )). 

Since the straightforward factorization cannot be used for baryons, we will apply 
the fourth sum rule to evaluate the operator Od in A;,. The fourth sum rule for it 
takes the form 

-^(A,|27ra,Oz,|A,) = -3 5] i^^ja^p . (89) 

The states appearing in the r.h.s. are orbital excitations of the A^ baryon with 
quantum numbers of the light degrees of freedom s^* = 1~. Their excitation energies 
are En and a'-"^ are the corresponding oscillator strengths describing semileptonic 
decays of Aj, to the analogous excitations of the Ac baryon; they are defined as in 



The first excited states appearing in this sum rule have been identified as the 
doublet of negative-parity baryons A+(2593) and A+(2625). Unfortunately no ex- 
perimental information is available to date on the transition amplitude a^^^ governing 
the decays of A;, into both of them, although it will be ultimately measured. 

The important piece of information would be the slope of the elastic IW 
function in A^. This quantity is more accessible than a*^^^ and will be measured in 
the near future at LEP. In the absence of the data we can use the second sum rule 
(Voloshin's "optical" sum rule) for A^^ = — rrih — Ma^ — rric. As discussed 
in |]33| (for earlier application see also |^2|, |38|), we can estimate pfj using just the 
excitation energy of the first states. We simply take Ma* ~ 2.615 GeV and the 
first excitation energy Ai ~ 330 MeV. Assuming A^ ~ 600 MeV and, therefore, 
Aa, ~ 900 MeV, we then have 

(p^^^ ^ ^A?Aa, ^ 0.15 GeV^ (90) 

(A similar estimate (/U^)aj, — | ^i^A,, — 0.45 GeV^ agrees well with the mass rela- 



tions 1 3^, for charm and beauty in the meson and baryon sectors). We note in 
passing that, most probably, the large mass of the heavy baryon implying larger A 
compared to B meson is due to larger slope p\ ; the higher- dimension operators, 
therefore, can be even smaller than in mesons. 

The sum of the expectation values for all three light flavors 2A' + A'^ is related to 



the expectation value of the Darwin operator [31|: 



Hence, we estimate the S'?7(3)-singlet color-octet expectation value in A^, as 

3 

47r 



K+,+s ^ |:(p?>)a ^ 0.036 GeV^' (91) 
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with uncertainty about 30-40%. The estimate of can be obtained if we neglect 
the small contribution of the non- valence strange quarks: 



A' 



A' ~ — 



0.018 GeV^ . 



(92) 



We note that we get a reasonable agreement with the quark model relation A' ~ A 
between the straight and octet expectation values in A;, for our central estimates, 
Eq. (|65|). It is interesting that the corresponding value of the diquark density at 
origin appears close to our central estimate for mesons —ooy (but larger than the 
alternative analogue of |\E'(0)p in mesons /^Mb/12). It also exceeds the estimates 
obtained in the QCD sum rules or quoted from bag models these analyses 
determined the combination |A' 
0.004 GeVl 



^A which generally emerged in the ball park of 



5 Nonfactorizable pieces in the matrix elements 
of the four-quark operators 

As was mentioned earlier, there are four operators (for a given light quark flavor) 
determining the corrections to the mesons widths. These are color-straight O and 
color-octet T. Each of these can contain either timelike (vector Oy, Ty) or spacelike 
(axial Oa, Ta) components of light and heavy quark currents. 

The color-octet expectation values r in general can be estimated using vacuum 
factorization, since the operators T coincide with the operators colorless in the s- 
channel up to 1/Nc terms (see Eq. (p3D). For such operators vacuum factorization 
is expected to work up to 1/Nc corrections. This gives 

^ JUpWb m _ _ ] ^ _o.015 GeV^ (93) 




TA = (^1-— J ^ 0.045GeV^ (94) 

It is interesting to note that the leading 1/Nc corrections to the factorization 
approximation can be estimated in a phenomenological approach. For this the ex- 
pectation value of the color singlet Osingi in (^) is written as 

{H,\{bTq){qTb)\Hh) = {H,\bTq\0){0\qTb\H,)+J2{H,\bTq\n){n\qTb\H,). (95) 

n 

The leading corrections to the vacuum factorization approximation are of order 1 
and come from one-particle intermediate states like vr [t], t]'), p (u), ai for B mesons, 
or light baryons for A5. The corresponding transition amplitudes have been evalu- 
ated in the QCD sum rules and lattice simulations with an accuracy sufficient 
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for determining the scale of the effects. Alternatively, they can be approximately 
obtained from the corresponding decays of charmed particles. 

Large- A^c arguments per se do not ensure that the vacuum factorization approx- 
imation works in the case of the color-straight operators, for the nonfactorizable 
contribution appears at the same order in Nc as the factorizable one. Their Wilson 
coefficients are not suppressed compared to those of the color-octet operators T (see 
Table 1), and they can be important even if their expectation values are formally 
subleading in l/N^. Moreover, the factorizable part of the expectation values has 
only a specific Lorentz structure which is subject to the strong chirality suppres- 
sion ~ m?^/m1 in the effects of weak annihilation (WA) in mesons. Nonfactorizable 



contributions there can be dominant ITT 



Nonfactorizable effects also appear as the expectation values of the non-constituent 
quark operators. Although they do not split the widths of and , they can 
differentiate the meson vs. baryon lifetimes. Numerically they seem to be strongly 
suppressed, with a possible exception of the Darwin operator which will be discussed 
below. 

The nonfactorizable effects in B mesons were first discussed in the framework of 
the 1/m expansion in |TI[ where the parametrization 



^(5|67^(l -75)gg7^(l -75)6|5) = ^^^^{v^v^v^ - g,g^^) (96) 



2Mb 



PbMb, 



^{B\h,{l-l,Yqqiu{l-l^)t%\B) = l^(voV^v, - goQ^,) (97) 

was suggested motivated by the analysis of the WA effects: neglecting the c quark 
mass WA is governed by Qo and Qs- (In the factorization approximation = 1 and 
Vo = go = Os = These parameters are related to u, r in the following way: 

UMbVs = -J-^v + - 4ry + ^TA (98) 

11 2 

flMBQs = -Jf^V - - 2ry - -TA (99) 

PbMbVo = -(l-^]u^V + l(l-^]u^A + ^ry - ^TA (100) 




3N. 



c 



flMsgo = -^^[^-irA^y - ^a[^-ij2)^A + j^ry + ^JA. (101) 

The inverse relations expressing uj and r via v and g are given in Appendix 3. 

The color counting rules suggest that Ty^A ~ while ujy^A ~ N^- The factor- 
ization estimates for Ty^A in the large- A'c limit are expected to hold with the 1/A^c 
accuracy. Therefore, knowledge of the color-straight operators allows to estimate 
the leading, 1/Nc terms in Vo and Qo'- 

Vo ^ --—^ — — — TTTT (valence) (102) 
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(The term —l/[2Nc) is absent for non- valence expectation values.) 

Since the non-valence expectation values appear to be suppressed, we only con- 
sider the valence matrix elements generically denoted by the superscript (f). Let 
us consider for definiteness the charged B meson; the corresponding parameters are 
then defined by Eqs. (p^),(p7D with q = u. Although uy and lja are not precisely 
evaluated, we still observe a clear tendency to cancellations in go and, therefore, 
suppression of the effects of WA. Say, for the exponential ansatz we get 

v^J^ ^ -0.07, g^J'^ ^ 0.004. (104) 

For the two-pole ansatz representing the upper limit in our estimates, we get 

vi""^ ^ 0.4, g^"^ ^ 0.05. (105) 

The nonfactorizable octet parameters seem to emerge suppressed. In particular, the 
expectation value of the operator responsible for WA is very small. The color-singlet 
expectation value gs was not estimated in the literature. It is natural to think W2 



that the scale of gs does not exceed that of go- The above estimates then illustrate 
the degree of suppression of the effects of WA when the c quark mass is neglected. 

It is appropriate to note at this point that there is a convincing experimental 
evidence that WA is indeed strongly suppressed in heavy mesons. The width dif- 
ference between Dg and D° is very sensitive to WA. Even though the literal l/rric 
expansion in charmed particles is hardly applicable at the quantitative level, the 
significance of such effects would have led to a large to^- tdo difference. Barring 



accidental cancellations one gets a typical estimate [|12 

\go, gs\ ^ 10"^ . 

We note, therefore, that our estimates, whatever tentative they are, indicate a strong 
enough suppression. It is interesting that the QCD sum rule estimates of the pa- 
rameters g made in 1992 by V. Braun ^ yielded close values, a few units xlO~^. 
Later evaluations gave Vo ^ 0.05, go ^ 0.05 |^, and Vo ^ 0.1, go ^ 0.03 ||4^; they 
were simplified in many aspects, though. While the expectation value of Vo generally 
emerges of the order of 0.05, our estimates for go seem to predict typically somewhat 
smaller values ~ 10~^, in a better agreement with the experimental indications. 

The non-valence expectation values are probably even further suppressed. Our 
estimates yielded w^"^) ^ (0.5 to 2) ■ 10"^ ^nd ^ -(0.25 to 1) ■ 10"^ with the 
dominant part coming from the axial current via the anomalous terms. 

For baryonic expectation values there is no vacuum factorization approximation. 
This does not mean, of course, that they are suppressed. The color-straight ex- 
pectation values A were estimated in the preceding sections to vary from 0.007 to 

^Private communication to N. Uraltsev in March 1992. Unpublished. 



31 



0.045 GeV^; the interval above 0.03 GeV^ seems improbable, though. The estimate 
of the color-octet ^ based on the evaluation of the Darwin operator yielded about 
0.018 GeVl 

A different parametrization of the valence expectation values was used in [E3|: 



2M, 



^ 'B-\ih,il-j,)u)iu^,il-j,)b)\B) = ^Ml^B^i^ (106) 



B 



^ -{B-mi-^,)u)iuil + ^,)b)\B-) = iMl^BM (107) 



2Mb /V V 2 

^ {B-\{h,il-l,ru)iu^,il-^,rb)\B-) = ^i^^e.if^) (lOJ 



2M, 



B 



^ {B-mi-^,ru){u{l + ^,)t%)\B-) = ^^^^j^^ e2(/i) . (109) 



2Mb ^ 
These parameters are related as follows: 

/|Mb B, = ~PbMb{v, - Ag,) = 4(ry + ta) + ^{uy + ua) (HO) 

flMB B2 = ~PbMb{vs - 9s) = ~2{ry - ta) - ^{ujv - ^a) (HI) 

PbMbSi = pMBivo-Ago) = -^{tv + ta) + (l " ^) {ooy + ooa) (112) 

/|Mb£2 = PbMb{vo-9o) = Y^i^v-TA) - ^(i-]^)k-^a). (113) 

The above estimates for the octet expectation values neglecting l/N^ terms look for 

ei^2 as 

£1 ^ -0.085 to 0.17 , 62 ^ -0.07 to 0.33, (114) 



while the QCD sum rule calculations read ei ~ —0.15, 62 — and ei ~ —0.04 it 



0.02,^2 ^ 0.06 ± 0.03 42 



For convenience, we give in Table 1 the central estimates of the four-fermion ex- 
pectation values discussed above in B mesons {ujv,a, 'Tv,a) and in (A, A' for a fixed 
flavor, u or d). Since the non- valence expectation values are strongly suppressed, 
we do not quote them here. 

Concluding this section, we note that there are two exact positivity constraints 
on the expectation values of the s-channel colorless operators: Vg—gs > 1 ivs—Qs > 
for non- valence) and Qs > 0. 

The first inequality follows from inserting a complete set of intermediate states 
\n) in the matrix element (67^(1 — 75)9) (^7°(1 — 75)^)1-8). We obtain 

= 1 + 1277^ E/ dfi{n)\{n\q^'{l-^,)b\B)\'>l, (115) 



B n 



32 





UJv 


lua 




Ta 


Au 


K 






Sects. 4.1-4.2 


-0.020 


0.045 






0.020 








factorization 


-0.011 


0.034 


-0.015 


0.045 






0.10 




4*^^ sum rule 






-0.028 






0.018 


0.18 


0.15 



Table 1: Estimated valence expectation values in B and A;,, in GeV^; factorized 
contributions assume = 160 MeV. 

where d/x(n) stands for the phase space. The 1/Nc contributions in the r.h.s. come 
when the intermediate states n are vr, p, ai, etc. For non- valence quarks the vacuum 
factorization contribution 1 in the r.h.s. explicitly showing \n) = |0) vanishes.^ In 
terms of the parameters Bi, Ei, the constraint ( |115D reads B2 > 1. Estimates of |^ 
give values for B^ compatible with 1. 

The second inequality is obtained by taking spacelike fi = u = i in (^). Sum- 
ming over i yields 

9s = r^i^E / df^in)Y.\{n\qYil-l,)b\B)\' > 0. (116) 

(For the B parameters this is -B2 > -Bi-) 

A similar inequality can be obtained for the baryonic matrix elements: 

A'-^A = ^X:/ dMr^)|(n|g(l + 75)&)|A,)r>0, (117) 

where we used the identity 

(A,|(6(l-75)g)(g(l + 75)&)|A6) = -^(Afe|0v|A5) - (A,|Tv^|A,) . (118) 

In the constituent quark model the bounds for B mesons become equalities; the 
relation Eq. (|117|) merely expresses the fact that the diquark wavefunction at origin 
is positive. It does not seem to be very restrictive for our estimates. There is an 
additional constraint on the expectation values of the operators with chirality flip 
for light quarks, however we are not interested in them. The above bounds can 
be refined and even the actual approximate estimates can be obtained evaluating 
the contributions of a few lowest intermediate states in the hadronic saturation 
Eqs. ([TTg)-([lTl. 

^These inequalities assume a physical regularization scheme for composite operators in which, 
for example, there is a power mixing of the four-fermion operators with the unit one, bb(0). For a 
recent discussion see, e.g. |3|, |35[ . 



33 



6 Corrections to the decay widths 



In this section we give the expressions for the corrections to the widths in terms of the 
effective four-fermion operators normahzed at a low scale. These expressions were 
originally derived in We present them here for completeness and book-keeping 
purposes, in a more convenient form. 

Let us introduce the notation AF for the operator describing the corrections to 
the inclusive width AF// of the beauty hadron Hb- 



AT 



2Mh, 



■{Hh\Ar\Ht) . 



;ii9) 



In what follows we neglect the effects generated in the KM suppressed decays which 
have a factor |Ktfe/KfeP, and by Penguin operators in T-C^eaki^B = 1) at the scale 
nib. 

Without accounting for the perturbative QCD effects in the domain k -C one 

has 



AT 



\V,b\\l - y? I U + cl + ^c,C2 ) m + Ol] + 4C1C2 [T^ + n] 



+ 2{2ciC2 + N^cl) 



C? + C2 + j;fCiC2 



{i + y)o^ + l{i-y)oi 



[I + y)T^ + 1(1 - y)Tiy^ 



Ol. + -{I - 4y)0% 



c? + ci + ^+ 



r^' + i(i-4y)ri' 



2Nr 



+ 2{2ciC2 + Ncci + ^ 



(120) 



rjic 

rpC I -'-A 



We denoted here y = ml /ml and d' = d cos 6^ + s sin 9c, s' = s cos 9c — d sin 9c. The 
TTic-dependence for the operators with external c quark legs is completely neglected. 
(Eventually they will lead only to Penguin-type operators which are estimated, ba- 
sically, in the leading- log approximation.) We included the contribution from the 
semileptonic decays with rii = 2 species of light leptons (the r contribution is sup- 
pressed by the phase space). Since numerically m^/ml ^ fJ'had/^b, keeping m^/mf 
corrections apparently is not legitimate in practice at all. We retain these terms 
only for getting an idea of the scale of the finite-mc corrections in the coefficient 
functions. 

The perturbative evolution below mj, in the LLA is particularly simple in this 
basis: the color-straight operators O do not renormalize. The color-octet operators 
T renormalize in a universal way with 77- = -iNc, except for the flavor singlet vector- 
like operator similar to Od which has anomalous dimension 7^, = 3Nc — |n/ where 
n/ is the number of open flavors {■jd = — yA'c). At the scale below the charm 
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mass the operators with the c-quark fields merely vanish. As a result, at the low 
normalization point /i we have 



(1-2/)' 



Ar 



01 + o^. 



Glml 

1 + 2/ 



1-2/ 
6 



oi 



2(1-2/) 



6(1-2/)' 



(l + y)r^' + 



+ C(l - 2/)' [4ciC2 (T^ + T;^) - (2ciC2 + N,cl) 
^/^^ , , (l-4j/)^ 



1-2/ 



where 



c 



2,flo-4/3 



(1-2/) 



^^(mc 



rjiS 

2^V 




3(1-2/) 



/?0 



rpS 

2^ A 



11 



-AT, 



9 

27rcD7— 
zvr 



9 



(121) 



'122) 



We wrote the contribution of the Darwin operator separately although it is related 
to the sum Ty + Ty + Ty. In this form cd emerges from the Penguin- type diagrams 
while the other terms do not include the annihilation diagrams. 

The coefficient cr, of the Darwin operator Eq. (BBT) takes the following LLA form: 



1 



Cd 



27ras(/i) 



4ciC2 - (2ciC2 + N^cl) [l+y 



C (i-y)^ 



(l-Z/)^ 




, ^(e-l) 



2ciC2 + N^cl + 



where 



2{n^+l) 

'a<j(mc)\ ^'^0-2 



V 



f 

■ W 



nf = 3. 



(123) 



(124) 



In principle, there is another source of the Darwin term in the width which comes 
from the l/rrif, expansion of the expectation value of bb and from the non-logarithmic 
terms in the expansion of the transition operator. They were calculated for the case 
of the semileptonic width in and We can estimate this correction to the 
LLA neglecting the deviation of the color factors ci , C2 from their bare values 1 and 
0, respectively, and neglecting the mass of the quarks (leptons) produced by the 
virtual W boson. In this approximation the possible effect from the ud (cs) loop 
cancels and we can use the calculations for the semileptonic widths. This yields to 
the leading order in «<, 



6c 



■D 



^^^^ (77 - 88y + 242/2 _ gy^ _ Sj/^ + ?>Qy^ In y) 



(125) 
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The overall non-log term in cd appears to be of the opposite sign to the LLA result 
and is roughly a half of it in magnitude. We conclude that the LLA estimate is 
accurate within a factor of 2 being, probably, on the upper side. We use the LLA 
expressions for numerical estimates below. 

At the next-to- leading order the Ai? = 1 weak decay Wilson coefficients ci(mf,), 
C2(m5) are | 

Ci = L13, C2 = -0.29, (126) 
corresponding to as{Mz) = 0.118. This gives 

2 

cl + cl + —ciC2 ~1.15, 2ciC2 + Nccl ^-0.40, 

4ciC2 ~ -1.3 , 2ciC2 + N^cj - 3.2 . (127) 

With this input the resulting values for the coefficients in Af are given in Table 2. 
To illustrate the uncertainty associated with the LLA we quote two sets of values 
corresponding to using as{mb) = 0.3 and to as{mb) = 0.2 (the former option can 
represent the choice of the ^-scheme strong coupling in the LLA expressions, which 
seems more appropriate to us). The coefficients A and B are defined as 

Af = (aIO^ + At.O'^ + Al^O^; + A\Ol + A\0i + A\Oi+ 

Z7l ^ 

B^T^ + B^T^ + B'yT{: + BlTl + B^Tf + B^T^ + 27rcz5^) (128) 

(i.e. using the 'redundant' basis including the Darwin operator to show explicitly 
the loop contributions). We also quote the values of i?y. By and By given by 

5^ = - 27ra,CB , B'^ = B'(,-2TTa,CD, B'y = B'y - 2nasCD ■ (129) 

The Cabibbo mixing is neglected here. 

For numerical estimates we use the values of the running low-scale masses rub — 
4.6 GeV, rric — 1.25 GeV and normalize the width correction ATh^ to the semilep- 
tonic width which is reliably evaluated in the OPE (for a review, see [^). The final 
estimates for these corrections then read as 



^£J!±^J_^£3l ^ 36 



^ 0.02 GeV^ 0.02 GeV^ 0.02 GeV^ 



0.02 GeV^ 0.02 GeV^ 0.02 GeV^ 0.02 GeV^ 0.02 GeV^ 0.02 GeV^ 

+ B^^ZI^ + S^^Zl) + BJ^ZI) + 0.8 {27rcD) , 
0.02 GeV^ 0.02 GeV^ ^0.02GeV^ 0.1 GeV^ 



^Note that these NLO values of ci.2 are immediately reproduced in the simple LLA if one uses 
the more physical TZ-scheme as coupling |471| . 
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asirrib) 


AU 


At 




Ai 




TDU 

BvA 


Bi 


B'v 


Bi 


B% 


0.3 


0.98 


-0.53 


-0.48 


-0.15 


-0.11 


-2.1 


0.70 


0.64 


0.20 


0.15 


0.2 


0.98 


-0.53 


-0.48 


-0.15 


-0.11 


-2.6 


0.84 


0.77 


0.24 


0.18 





27rc/3 


B^ 




B'v 


0.3 


-0.80 


-1.3 


1.5 


1.4 


0.2 


-0.56 


-2.0 


1.4 


1.3 



Table 2: Values of Wilson coefficients for the total width 



0.36 



.02GeV' .02GeV' .02GeV' '".02GeV^ '".02GcV^ "'.02GeV^ 



^.02GeV^ ^.02GeV^ .02GeV^ ^.02GeV^ .02GeV^ .02GeV^ 

(130) 

where {Oy) = ^tt — {Hb\Oy\Hi,) , etc. We recall that the expectation values in B are 



denoted by u for color-straight operators O and by r for the octet ones T, Eq. {paj)] 
for Ab these are A and — |A', Eq. (^). 

It is interesting to note that, regarding the Nc counting rules one can view the 
Wilson coefficients of the color-straight operators to be while the coefficients of 
the octet operators as l/N^. This is true if we recall that formally Ci(mb) = 0{1) 
while C2(mft) = 0{1/Nc). These are not mandatory assumptions for the large- 
Nc analysis: smallness of a particular perturbative renormalization can always be 
compensated by large logarithms of Mvy/m?,; in any case the nonleptonic weak decay 
coefficients Ci^2 are external to QCD itself and can be taken completely arbitrary. 
Nevertheless, their numerical values fit well such a naive assignment. 

Our procedure of evaluating the l/ml corrections to the widths then gets justi- 
fication in the formal Nc counting rules: we take at face value the color-straight 
expectation values appearing with the coefficients ~ Nj^, and take only the leading 
factorizable values ~ Nc for the color-octet operators which come with the sublead- 
ing coefficient l/Nc- This formally sums all leading corrections ~ N^ in the decay 
widths. 



7 Discussion 

We have considered the expectation values of the four-fermion operators which are 
encountered in the l/mq expansion of the inclusive widths of beauty hadrons. The 
size of the color-straight operators used to be most uncertain in B mesons, since the 
factorization approximation a priori is not expected to be accurate for them. On the 
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other hand, just these operators have the most direct meaning being analogues of 
the usual wavefunction density |\f (0) p. Using the exact relation of their expectation 
values to the momentum integral of the elastic transition amplitudes, we estimated 
these expectation values employing reasonable assumptions about the behavior of 
the formfactors. We showed that the actual large-g^ asymptotics of the light quark 
amplitudes in heavy hadrons is l/{q^Y^'^ rather than as has been believed 
based on simple-minded quark counting rules. We also calculated the anomalous 
dimension of the color-straight operators and their mixing with the octet operators, 
the effects absent at order ag- The order-cc^ corrections appeared to be numerically 
enhanced. 

In our estimates of the valence expectation values their size obtained from the 
two-pole ansatz can be considered as an upper bound. A more reasonable expo- 
nential approximation which suppresses the contributions of momenta above 1 GeV, 
yields smaller results. We accept it as a typical lower bound for the color-straight 
expectation values. Although the accuracy of the central estimates cannot be too 
good, they probably hold better than within a factor of two. 

Our estimates, in principle, include a source for non-valence expectation values. 
It is related to a different g^-behavior of formfactors describing different isospin 
amplitudes at < 0. We have it mainly as the different masses of the isosinglet 
resonances saturating the formfactors in the t-channel, compared to the correspond- 
ing flavor nonsinglet particles {i.e., annihilation shift of masses). Except for r/', 
experimentally these splittings are rather small, and our literal estimates thus yield 
a strong suppression. We are not sure if this really applies to the color-straight 
expectation values; the actual suppression can be softer. 

We observe a weaker suppression of the non-valence color-straight matrix ele- 
ments for the operators with the axial current. It is related to the nonperturbative 
'annihilation' effect, in particular, the axial anomaly in QCD and its solution of the 
U{1) problem. We conjecture that the dominant effect is the mass shift of the lowest 
pseudoscalar state rj' while the splitting of the massive resonances (in particular, ax- 
ial) or the effect of the possible difference in the singlet and triplet couplings (ja(0) 
and G^'* (0) is smaller. Then we get a tentative relation 

^(fi1^„...J.„..<^.h„7.<..w|BT . (131) 

This estimate has the correct scaling Numerically, the axial non- valence ex- 

pectation values appear to be suppressed by a factor about 0.1. We note that the 
numerical suppressions of various non-valence effects typically is stronger than the 
naive factor 1/3 which can be expected if their justification is merely the large 
A^c = 3. 

An interesting indication from our estimates is that the possible nonperturbative 
vitiation of the chirality suppression of WA in B mesons emerges at a rather low level 
(it is governed by the combinations iujv + \'^a) , [tv + |ta) ) ■ For the color-straight 
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operators (where the effect a priori can be significant), the hteral suppression is by 
more than an order of magnitude, in accord with the evidences from charmed mesons. 
In our approach the origin of the suppression roots to the fact that —Ga{0) ^ 
1/3. The WA effect of the octet operators can be probed in the difference of the 
semileptonic b —>■ u distributions in and JlT] . 



The chirality suppression of WA can be eliminated already in the perturbative 
evolution of the effective operators. This does not happen in the LLA [^]. Our 
NLO calculations show that it does not happen at this level as well. It is interesting 
to check this property for the two-loop diagonal renormalization of the color-octet 
operators. In any case, we expect it to be lifted in three loops; also, the non- 
logarithmic gluon corrections at A; ~ rrib defining the initial values of the Wilson 
coefficients must generate the chirality non-suppressed effect at some level. 

Let us now turn to the phenomenological consequences of our analysis. The 
estimated expectation values are typically of the order of, or somewhat larger than 
the factorization values (when the latter are possible) at Jb = 160 MeV (the fac- 
torization value of LJv is additionally suppressed, and our estimates only partially 
reproduce this). The actual expectation values of the color-straight operators can 
be smaller if, for example, the formfactors change sign at — <^ 1 GeV^. Such sub- 
tleties are not properly captured by the simple models we relied upon. On the other 
hand, larger values than quoted in Table 1 are improbable, at least if the nonper- 
turbative dynamics we account for are dominated by the momenta not exceeding 
IGeV. 

The relevance of the latter assumption for the analysis of the inclusive widths is 
easy to see, say, on the example of the effect of interference (dominant in B mesons). 
The decay rate of the process b ^ uj: + (cd)^ is proportional to = {pb — kY. 
At fc^ = one has = ml — 2pi,k, and this constitutes only about 12 GeV^ vs. 
ml ~ 21 GeV^ already for |fc| = 1 GeV. At the same time the usual relation of the 
1 / ml effects via the expectation values of the corresponding four-fermion operators 
assumes that = mf. Therefore, if \k\ becomes as large as taken above, the 
validity of the leading-order expressions breaks down. In any case, accounting for 
the effects like interference in the usual way is legitimate only if their impact is much 
smaller than the partonic width of a particular quark channel. It is worth noting, on 
the other hand, that the assumption that the nonperturbative contributions to the 
expectation values come from momenta not exceeding 1 GeV is built in the approach 
of the QCD sum rules. 

At first sight, WA in mesons and 'weak scattering' (WS) in baryons can get en- 
hanced, in contrast to interference, if the quark momenta saturating the expectation 
values of the operators are large. Such a conclusion, even though eventually may 
prove to be correct, cannot be justified a priori, and even the sign of the correspond- 
ing corrections to the standard expressions is not known. All such effects manifestly 
go beyond the 1/m expansion truncated after l/mf terms. For this reason, simply 
assuming large expectation values in B particles does not allow one to boost signif- 
icantly the lifetime differences respecting the self-consistency of the simplest 1/m;, 
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expansion. 

Bearing in mind all reservations made above, we still quote the central values 
for the corrections to the inclusive widths stemming from our analysis: 

- 0.36(-l.lpi - 1.2d) , ^ ~ 0.36(-0.15wA - 1.2d) , 

J- si J- si 

^ ~ 0.36 (2.2ws - Ipi - Id) ■ (132) 

J- si 

Here we showed separately the effects of different light flavors: of the operators 
{hh){uu) responsible for PI in B and WS in A^, and of {bb){dd) generating WA in B 
and PI in Af,. We singled out the contribution of the Darwin term. Even though it 
may seem to be a computational separation, it is legitimate, for it can be formally 
carried through the dependence on the number of light flavors. Being a flavor singlet, 
the Darwin operator does not differentiate the lifetimes of charged and neutral B 
(also of Bs to the extent that S'f/(3)fl is a good symmetry). 

The above estimates generally support the original theoretically predicted pat- 
tern of the lifetimes. The non-valence effects seem to be strongly suppressed. The 
main effect is destructive PI in B~, about —4%, while WA is small, at a half per- 
cent level. Moreover, literally we get the effect of WA decreasing the width, the 



possibility originally discussed in [^8|, |TT|, and which may seem to contradict 
the naive interpretation of WA. The overall difference of r{B~) and r{B^) appears 
about —4%. The major effect is WS in A;,, 8.5%, but it is partially offset by in- 
terference, —3.5%. The difference between r(Aft) and r{B^) is literally 6%. These 
estimates fall close to the expectations quoted in the review We note that the 
often discarded Darwin term (e.g., in [^]) typically decreases the width by about 
4%, although literally we get its effect in B and A^ close to each other. Including it, 
the overall decrease in the A;, lifetime from the four-fermion operators at the order 
1/ml comes out only at a percent level while tbo increases by 5% and tb- by 9%. 
The overall absolute shift is not too interesting by itself though, since it depends on 
the exact definition of the parton width. 

It is worth noting that the corrections we addressed do not formally exhaust the 
1 / ml terms in the asymptotic expansion of P^^ - they come implicitly as well from 
the expectation values of the kinetic and chromomagnetic operators which appear 
at the level of corrections. These expectation values in the actual b hadrons 

differ from their asymptotic values at oo by terms ~ 1/m;, P]. In particular, 

these deviations contain the expectation value p^g of one new local heavy quark 
operator, the convection current (or spin-orbital) one. (This operator cannot appear 
independently in the expansion of the transition operator describing the inclusive 
width since it is not Lorentz-invariant.) These corrections do not affect P^- — Pbo 
but, in principle, are present in P^ — Pa^. Their practical neglection nevertheless is 
legitimate: such effects are included in the existing uncertainty of the differences of 
the expectation values /i^ and /i^ of the D = 5 operators between B and Af,. So far 
these expectation values are estimated without considering corrections to the heavy 
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quark limit; for example, the value of fiQ in A;, is nonzero but generally of the order 
of AQf^j^/mfo. All such effects are also expected to be numerically insignificant. Let 
us recall that in B mesons the pj^g expectation value is suppressed to the extent 
that their two-particle description is applicable 0. 

Our analysis does not indicate a crucial impact of the nonfactorizable contribu- 
tions in the low-scale expectation values on the B lifetimes conjectured in |^3| or 
later speculations that r(i?+) can even exceed r{B^) by a significant amount. 

The small experimental lifetime of A5 thus remains a challenge for the straight- 
forward l/mj, expansion. An accurate measurement of the semileptonic width of A^ 
(or BRsi(Aft)) would help to shed light on the origin of the problem. The gap be- 
tween the experimental value of ta^ and the theoretical expectations could have been 
reduced by a significant enhancement of WS and suppression of PI in A^, according 
to the natural guess about the role of the spectator momentum we mentioned above. 
Since these effects originate from the quark decay mode b — > cud constituting about 
60% of the total width, a 25% effect in the lifetime would signal a more than 50% 
enhancement of this channel. Clearly, such an effect is not possible for the spectator 
quark occupying only a small fraction of the total phase space in the decay, and 
would require non-conventional composition of the heavy hadron. The standard 
calculation of the 1/ml terms neglecting the effect of finite spectator momenta is 
not applicable for quantitative description of such large corrections. For example, 
the expectation values of the Darwin operator would be in general much larger, 
likewise the mass scale governing the size of higher- dimension operators for 
and higher-order corrections must be higher in this situation. 

Note added: When this paper was prepared for publication, a new improved 
QCD sum rule calculation of the four-fermion expectation values appeared pU[ ; 
the quoted results correspond to Vo — —0.03, go — 0.003, and Bi = 0.60 ± 0.01, 
B2 = 0.61 ± 0.01. It can be suspected, however, that the stated small errors did 
not adequately reflect the uncertainties inherent in the determination from the sum 
rules per se. 
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Al Combinatorial relations 



Here we quote two general algebraic relations which are useful in calculating the 
renormalization of amplitudes containing static heavy quarks. 
For any set of N numbers Xi, ... ,Xn 



k=o \ 2^1=1 
(it is assumed that H^^ = 1 if m > ^2), and 



fc=0 \ z^«=i -^1/ \ l^l=j -^i , 




(Al.l) 




The proof will be given below. 



. (A1.2) 



The sums of the type ( |ALl| ) are reminiscent to those appearing in calculating 
the renormalization of any color-straight operator of the type bb ■ Ought- The sums 
similar to Eq. (|A1.2|) emerge in calculations of mixing of an arbitrary heavy quark 
operator into the color-straight operators, bTb ■ Ought bb ■ Ought where T is any 
color matrix. 

For the color-straight weak vertex bb the product of color matrices on the heavy 
quark line does not depend on the location of the weak vertex in respect to the 
gluon vertices. The fc-th term in the sum Eq. ( |A1.1| ) corresponds to the diagram 
where the first k gluons attach to the initial b quark while the last N — k gluons 
attach to the final-state quark, Figs. 4a,b. We thus do not sum over permutations 
of gluons (their time ordering is fixed) but combine + 1 possibilities to place the 
weak vertex. The analogues of are Uk, the energies of gluons flowing into the 
quark line. With this identification the structure of the product of the heavy quark 
propagators is reproduced. 

In dressing a non-straight operator the gluon and weak vertices do not commute 
and moving the weak vertex would change the product of color matrices. However, 
calculating mixing into the color-straight operators amounts to taking trace over 
color indices of the initial and final state quarks. Then one, instead, can perform 
a cyclic move of the leftmost gluon in the initial state to the latest position in the 
final state, and vice versa. Figs. 4c,d. In considering Eq. ( [A1.2D we thus imply 
combining all graphs obtained by the cyclic permutations of a particular diagram. 
Both considerations apply for any color representation of the quarks and gluons. 

Taken naively, the relation ( |A1.1|) would suggest that the renormalization of 



the color-straight operators vanishes to all orders (already in the sum of the above 
groups of diagrams, before actual integration over all gluon momenta). Likewise the 
identity (|A1.2|) would look like the property that the octet operators never mix with 



the straight operators. This is not so, however. The reason is that the identities 
Eq. ( |A1.1D and (|A1.2|) apply only if the external b quarks are exactly on shell so 



that their nonrelativistic energy vanishes, E = 0. In this case the quark propagators 
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J3 ^^'4 ^5 '-^2 ^4 ^5 '^l 

2 d 



Figure 4: The diagrams combined for the color-straight operators (a,b) and for 
the mixing into the color-straight operators (c,d). The solid box denotes a color- 
straight operator, the blob in the diagrams c and d stands for an arbitrary heavy 
quark operator. Only two of six (A^ = 5) diagrams to be combined are shown in 
both cases. 

generally become IR singular when integrated over the gluon momenta, and must 
be regularized by a small imaginary part —ie in each heavy quark denominator. 
Alternatively, this is done by a shift of the external heavy quark energies by an 
infinitesimal imaginary amount. This regularization translates into the shift of all 
xi which, however, is of the opposite sign for the initial-state and final-state gluons. 
For example, for the sum in Eq. (|A1.1|) 



xi ^ xi + ie ioT I < k and xi ^ xi — ie for I > k 



(and the opposite shift in the sum in Eq. ([A1.2|) ). This infinitesimal shift of de 



nominators leads to the fact that the sum of all diagrams does not vanish exactly 
but contains certain (5-functions of combination of energies corresponding to a cer- 
tain on-shell heavy quark inside the diagrams. Nevertheless this kills some of the 
integrations over u and simplifies the remaining integrals. 

Let us prove identities ([Al.l|) and (|A1.2|) . This can be done most simply by 



using the following trick. We can consider the sum as a rational function of the 
variable xn (for example), at xi, ... ,xn-i arbitrary but fixed. If we show that the 
residue of this function at any potential pole vanishes, this would mean that the 
whole function vanishes identically. 

For the sum in Eq. ( [Al.lD this is particularly simple. Presence of a pole means 
that at certain k some of the denominators with j = jo vanish, with either jo < k 
(to the left of the weak vertex) or jo > k (to the right of it). Let jo < k, for 
example, and therefore J2iLi = 0. Then the same vanishing denominator will be 
present for all diagrams corresponding to A; > jo, and it will change only for k < Jq. 
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Moreover, all terms with k > jo will have the common factor 



T-rio-l 



which is the product of the propagators to the left of the one which vanishes. 

The remaining factors will be different, but for k = jo+l? ■■■ ,N their sum exactly 
reproduces the l.h.s. of Eq. (|A1.1|) for the set of Xj^+i, ... ,xn (that is, the case of 

— jo gluons) owing to the on-shellness of the jo-th propagator (the condition 
X^ili ~xi = 0). The induction from the obvious case N = 1 immediately proves 
Eq. ( Al.l ) for arbitrary A^. 





^5 ^ ^4 



Figure 5: Graphic illustration for the case N = 8. 



The proof of the identity Eq. ( |A1.2| ) is a little more complicated. To phrase it, 
it is convenient to close the heavy quark line and map it onto the circle. Figs. 5. 
The weak vertex can be referred to as North Pole whereas the infinity can be called 
(with some reservations) South Pole. Every arc on the circle can be attributed the 
corresponding energy denominator. Proceeding from the k-th arc to the k + 1-th 
arc clockwise decreases the denominator by Xk. The values of all denominators are 
fixed by the condition that the arc containing the South Pole (the Infinity arc) 
has vanishing denominator (correspondingly, it is excluded from the product of 
propag ators in Eq. ([AT^ ). 

With this image it is easy to establish the vanishing of the residues in the sum 
Eq. (|A1.2|) as well. A pole would appear due to the vanishing of the denominator 
of some other arc with j = jo; it is indicated by the star in Fig. 5a (the Zero arc). 
This figure shows the case of jo = 6 and k = 4. It is easy to see that the residue 
is exactly canceled by the configuration with jo ^ k, that is, when the Infinity arc 
and the Zero arc are interchanged. Fig. 5b. 

Indeed, due to vanishing of the denominators at the both arcs all other denomi- 
nators in Fig. 5b are equal to the corresponding denominators in Fig. 5a. To get the 
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residue one must merely remove the two vanishing propagators from the product 
and take it with the factor —1 when the Zero arc is clockwise from the Infinity arc 
and with the factor +1 otherwise. This cancellation in Eq. (|A1.2 ) reads as 



JO 

E 

yl=k+l 



2^i=k+i Xl 2^i=j Xl 



jo 

E 

^l=k+l 



-Xl n 



N 

i=io+i 



JO 



■ n-^° - 



at 



l=jo+l " S«=j l=k+l 



-Xl 



Thus, both identities (|A1.1|) and ( |A1.2| ) are proved. 



A2 Two-loop anomalous dimensions 

For the order-a^ hybrid renormalization of the heavy quark operators QQ qq the 
identities discussed in Appendix 1 say that summing over all attachments of the 
gluon to the heavy quark line results in 5{uj). Therefore, the integration over d^fc 
cannot produce an UV logarithm since it would require an odd power of k in the 
integrand. This is not possible in the simple one-loop diagram. The one loop 
renormalization of the straight operators coincides, therefore, with that of the light 
quark bilinear, while the octet-to-straight mixing is absent. For the vector or axial 
currents we consider, the overall one-loop renormalization vanishes. For the octet 
operator an additional contribution to the diagonal renormalization comes from the 
gluon exchange between the heavy and light lines. 

In order a1 both the renormalization of the color-straight operators and the 
straight-octet mixing occur. We do not consider the 0{a1) diagonal anomalous 
dimension of the octet operators. Since the 0{as) one does not vanish, the 0{a1) 
anomalous dimension is scheme-dependent. For the light-quark currents we are 
interested in, only nonfactorizable diagrams must be considered where at least one 
gluon connects the heavy quark line with the light part of the diagram. 

The hybrid anomalous dimensions are given by a (single) logarithmic UV diver- 
gence of the diagrams in the limit mq — >■ cxd, \k\ <C mg. In the Feynman gauge 
we adopt for computations, only 18 "double exchange" diagrams where two gluons 
connect light quark line with the heavy quark line each, yield the log. All other 
diagrams where there is only one gluon vertex either on the heavy quark or on the 
light quark lines, are finite for symmetry reasons similar to the one-loop case, or (in 
the case of dressing the octet operator) yield only the octet structure we are not 
interested in. 

Combining the diagrams into the groups of three according to the rules described 
in Appendix 1 (all locations of the weak vertex on the heavy quark line for the 
color-straight operators, or cyclic permutations of the QQg vertices for the octet 
operators) we get, at fixed values of the gluon momenta fci, k2 the sum of the heavy 
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quark propagators in the form 

1 1 
— 2niS(uJi + UJ2) or — 2niS(uJi + UJ2) (A2.1) 

for the color-straight operators, or 

- 27riS(uJi + UJ2) — ^ — + 27ri (v—SM - V—Sm) (A2.2) 
cui + ie V cui UJ2 J 

(and uji <-> UJ2) for the octet operators. In view of the uj — > —uj symmetry of the 
integration only the structure —2t:'^8{ijJi)5{uj2) survives, and the resulting integrals 
contain simple purely three-dimensional expressions given below . By dimensional 
counting they all are logarithmic; they do not vanish since integrations runs over 
two spacelike vectors. 

Dressing of color-straight operators QQ qVq 

The six groups of three diagrams in turn fall into three types which differ by the 
location of the gluon vertices on the light quark line, Figs. 6a-c. Each diagram can 
have gluon lines twisted or not. Their expressions are 



h = yC" [r7^7o7^7o] j 



(27r)3 (2^)3 k^k2(j:^ + k2y 
h = y C [707.r7.70] j (^^4^ (A2.3) 

h = y C [7o7.7o7.r] j ^^li^^l^ • 
The color factors C are 

C, = [tH%[tH\ = 1 (^1 - [!],[!], - i- n[t\ (A2.4) 

for "twisted" and "non-twisted" diagrams, respectively. 

For r = 7o or 7075, twisted or non-twisted separately, we have 

4+Wc--^.cry^^^^^^-^-,,cr^y -. (A2.6) 

For r = 7i or 7^75 

r ,r ,r _l Arr f ^'^^ ^'^^ " ^^)' + ^^^^ " ^^^^^^ + + ^^fkl 



46 



Figure 6: Diagrams showing the different attachments of gluons to the hght quark 
hne. 



The sum of all diagrams for arbitrary F takes the form 

QQqTq ^ (1 + ^(1-2^) InA ) QQ qVq + f 1 - 4t 1 InAQfQ qtTq . 



(A2.8) 



Mixing of octet operators Qt^Q qt^-Vq into color-straight operators 

Taking the trace over the heavy quark color indices we likewise can combine the 18 
diagrams into 6 groups belonging again to the pairs, where each pair has the same 
location of the qqg vertices but different trace of color matrices along the heavy line. 
For example, for F = 70 the projection onto the straight operator yields 



{la + h + Ic) [straight - "^^ (^^3 + C,)TJ ?4?4.^ , C ^2 



= .^(C3 + C.)r^f f , (A2.9) 

where the color factors 

The same renormalization emerges for other Lorentz structures F as well. 

For the flavor-singlet operators additional, annihilation diagrams are possible 
where the qq line forms a closed loop. It is easy to see that for the vector current it 
does not contribute. If the operator is color- straight, only two gluons can come out 
of the quark loop. The analogue of the Furry theorem leads to the cancellation of 
the two possible diagrams. 
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For any color-octet operator the sum of the diagrams where one of the gluons 
connects the external light and heavy quark lines yields only the octet operator in 
analogy with the one-loop diagrams. All other diagrams can obviously produce only 
the octet operators as well. 

For the color-straight operator with the axial light quark current, both gluons 
must come out of the quark loop. The expression for the triangle subgraph does 
not differ from the Abelian case |51|. The sum of the diagrams where one of the 
gluons is attached to the external light quark and another ends on the heavy quark 
line yields only non-logarithmic contribution. The diagrams when both gluons are 
absorbed by the light quark legs describe the two-loop renormalization of the singlet 
axial current and differ from the classic Abelian result only by the color factor 

Using Eqs. ( |A2.8| - |A2l!0|) and the definition Eq. (|28D , we arrive at the 0{al) 
matrix of the anomalous dimensions given in Eq. (E^) and in the text following it. 



A3 Relations between parametrizations 

Here we collect the relations between different parametrizations of the expectation 
values of the four-fermion operators in B mesons. 



Hadronic parameters suggested in Ref. [|T^] are given by 

PbMbVs = 



IbMbVo 
flMB9o 



N, 



■u^v - ^u^A - 2ry - 



SAT, 



(A3.1) 
(A3.2) 

TA (A3.3) 



■\ f 1 - 4 - ^1 - 4 I + ^ry + ^TA . (A3.4) 



6 



3N, 



The inverse relations read as 

1 



(jJy 
uja 
Ty 

TA 



PbMb 



JbMb 
PbMb 



PbMb 



-l^o -9o- 

-Vo - Sgo + 



■Vo + 



AN, 



1 


1 ■ 


f s - 

AN, 


2nM 


3 


3 ■ 




2N,^'\ 






-J' 8 I 




3 3 3 / 1 \ 3 / 1 \ 



AN 



2Nr 



A 



(A3.5) 
(A3.6) 

(A3.7) 
(A3.8) 
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We recall that for valence quarks v^^'^^ = 1 while v^^^'^^ = gl^'^^ = g^^^^ = . 



For parametrization of 

P^MbB^ = P^Msiv, - igs) 
PbMbB, = pMsiv^-g,) 
PbMbE, = pMB{v,-4g,) 

PbMbS^ = pBMB{vo-go) 

with the inverse relation 

1 



4(ry + ta) + —{ujv + CO a) 
-2(ry - Ta) - -^{uv - uja 



-Wirv + TA) 



1 - 



(A3.9) 
(A3. 10) 
(cuv + uja) (A3.11) 



— iTv-TA)--{l- — ]icov-coA), (A3.12) 



(jJV 

loa 
rv 
Ta 



PbMb 
PbMb 



AN, 



1 1 

-D2 + -ei - £2 



c . 



2N, 



5i + ——52 + -ei + ea 



2Nr 



IIMb 



1 - 



iV2 



Bi 



1 



1 - 



B, 



B 



-1-— fil + - l-— fi; 



4 



AN, 

1 

AN, 



1 1 



■ei 



2N, 
1 

2N, 



■£2 



(A3. 13) 
(A3.14) 

(A3.15) 
.(A3.16) 



All these relations hold for each light quark flavor separately. 

In the Ai? = 2 transitions -8°^-, — > B^^^ determining the width splitting in the 
B-B systems one encounters two four-fermion operators 0, both color-nonsinglet 
in the t-channel. They are naturally parametrized as 



(S,|6,7^(l-75)g* 6,7.(1 -75)g^|5,) = -2 Pb [v P^P, - g g^^Mf,) (A3.17) 
The non- valence matrix elements vanish. There is a standard notation Bb for 



{Bg\baa{i - 75)g* 6j-7a(i - i5)q^B,^) 



-^i'^ + Jf] BbPbMI. (A3.18) 



The anomalous dimension of this operator equals two anomalous dimensions of the 
bq currents, so that -B^ is renorm-invariant in one loop 0, |^ (all operators above 
are normalized at the low point, not at rrih). The combination of the operators 
corresponding to the v structure also renormalizes multiplicatively in one loop; its 
anomalous dimension was calculated in 0. Power mixing of these operators is 
absent. 
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